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Abstract 



The squared Bessel process is a l-dimensional diffusion process related to the 
squared norm of a higher dimensional Brownian motion. We study a model of n 
non-intersecting squared Bessel paths, with all paths starting at the same point 
a > at time t = and ending at the same point 6 > at time t = 1. Our interest 
lies in the critical regime ab — 1/4, for which the paths are tangent to the hard edge 
at the origin at a critical time t* G (0, 1). The critical behavior of the paths for 
n oo is studied in a scaling limit with time t — t* + 0(n~^/'^) and temperature 
T = l + 0{n~-^^^). This leads to a critical correlation kernel that is defined via a new 
Riemann-Hilbert problem of size 4x4. The Riemann-Hilbert problem gives rise to a 
new Lax pair representation for the Hastings-McLeod solution to the inhomogeneous 
Painleve II equation q"{x) — xq{x) + 2q^{x) — v, where i/ = a + 1/2 with a > — 1 
the parameter of the squared Bessel process. These results extend our recent work 
with Kuijlaars and Zhang [T^ for the homogeneous case = 0. 

Keywords: Squared Bessel process, non-intersecting paths, Painleve II equa- 
tion, determinantal point process, Riemann-Hilbert problem, Deift-Zhou steepest 
descent analysis, correlation kernel (Christoffel-Darboux kernel), multiple orthogo- 
nal polynomials, modified Bessel function. 
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1 Introduction 

The motivation of this paper is the recent surge of interest in a model of non-intersecting 
Brownian motions at a tacnode. The model is illustrated in the third picture of Figure [U 
where we have two groups of Brownian motions which are asymptotically supported 
inside two touching ellipses in the time-space plane. The interest lies in the microscopic 
behavior of the paths near the touching point of these two ellipses, i.e., near the tacnode. 

The tacnode model was recently studied via different methods by several groups of 
authors. The model is studied in a discrete symmetric setting by Adler, Ferrari and 
Van Moerbeke [2j. A different approach in this case is due to Johansson [23]. Further 
developments are the study of a double Aztec diamond model by Adler, Johansson and 
van Moerbeke [3], and the non-symmetric tacnode by Ferrari and Veto [15]. Another 
model with a tacnode but with very different properties is discussed in [4]. 

In a joint work with Kuijlaars and Zhang [13] we study the continuous non-symmetric 
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Figure 1: n = 20 non-intersecting Brownian motions at temperature T = 1 with two 
starting and two ending positions in the case of (a) large, (b) small, and (c) critical 
separation between the endpoints. The horizontal axis stands for the time, t G [0,1], 
and the vertical axis shows the positions of the Brownian motions at time t. For n — )• oo 
the Brownian motions fill a prescribed region in the time-space plane which is bounded 
by the boldface lines in the figure. In the case (c), the limiting support consists of two 
touching ellipses which touch each other at a critical point which is a tacnode. 

version of the tacnode model. Our results in [T3] express the critical correlation kernel in 
terms of a certain Riemann-Hilbert problem (RH problem) of size 4x4, related to a new 
Lax pair representation for the Hastings-McLeod solution to the Painleve II equation. 
Recall that the standard RH problem for the Painleve II equation has only size 2x2 
|16j : see also |25j for a Lax pair with 3x3 matrices. 

Recently Duits and Geudens [H] use the 4x4 RH problem from [13j to describe a 
new critical phenomenon in the two-matrix model. Interestingly the kernels in [131 [T3] 
are built from (basically) the same 4x4 RH problem but in an essentially different way. 

The goal of this paper is to extend the approach of [13j to a hard-edge situation. 
To this end we consider a model of non- intersecting squared Bessel paths [1^, see also 
|3H [32| [271 [28t [29] . In the critical case we can get a situation where the limiting hull of 
the paths is touching the hard edge at a certain critical time t*. This is shown in the 
third and fourth picture of Figure [2j We will refer to the touching point as a hard-edge 
tacnode. (This is an analogue of the higher-order tacnode in the classification scheme in 
|33[ p. 72].) We will see that the microscopic behavior of the paths near the hard-edge 
tacnode is again expressed by a limiting kernel defined in terms of a 4 x 4 RH problem, 
which is now related to the Hastings-McLeod solution to the inhomogeneous Painleve II 
equation. 
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Figure 2: n = 50 non-intersecting squared Bessel paths at temperature T = 1 starting 
in a > and ending in 6 > with (a) a = 6 = 2, (b) a = 6 = 0.25, and (c) a = b = 0.5, 
illustrating the case of large, small and critical ending positions, respectively. In the case 
(c), the paths are tangent to the hard edge at time t* = 1/2 (hard-edge tacnode) with a 
fourth order contact. A magnified picture near the hard-edge tacnode is shown in (d). 

Our scaling limit near the hard-edge tacnode will be more general than the one in 
[13] . In the latter paper we used a 0(n~^/^) scaling of the endpoints, or equivalently of 
the temperature T. In the present paper we will use an extra 0{n~^^^) scaling of the 
time t, similarly to [21 [T^ I24j . At the level of the RH problem this leads to an extra 
parameter r, as in [2]. While the extension of the RH problem with the parameter r 
was already known to the authors of |13j . the solvability of the extended RH problem is 
a quite non-trivial fact. This solvability was recently established by Duits and Geudens 
in the symmetric case |14j . We will make use of their ideas. We note that the solvability 
of the extended RH problem in the general non-symmetric case has not been settled yet. 

Now we describe our model in more detail, following [12j. Let X{t) be a squared 
Bessel process with parameter a > —1, i.e., a diffusion process (a strong Markov process 
on R with continuous sample paths) with transition probability density [26^ [30] 




x,y > 



(1.1) 



V. 



y>o, 



(1.2) 



(2r)°+ir(a + 1) 



e 2r 
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where r stands for the time, x,y for the position and where 



is the modified Bessel function of the first kind of order a; see jH p. 375]. If d = 2(a + 1) 
is an integer, the squared Bessel process behaves hke the square of the distance to the 
origin of a d-dimensional Brownian motion. Some apphcations of this diffusion process 
in mathematical finance and other fields can be found in [26] , see also [191 [28} [29t [M] . 

We will consider n independent squared Bessel processes Xj{T), j = 1, . . . ,n. It is 
then convenient to perform a rescaling of the time r i— )• t by 

r = Lt. (1.4) 

The proportionality parameter T > can be interpreted as the temperature. It has 
a natural physical interpretation, in the sense that high temperature corresponds to a 
large variance of the squared Bessel process and vice versa. 

We consider a model of n independent squared Bessel processes, starting in a > at 
time t = 0, ending in 6 > at time t = 1, and conditioned (in the sense of Doob) not 
to intersect in the whole time interval t £ (0, 1). This model was studied by Kuijlaars, 
Martmez-Finkelshtein and Wielonsky |3H [32] for 6 = and for general a,b > in 
the non-critical case in [12]. We also mention the large deviation principle recently 
established for a closely related model in [20] , 

Let us fix the temperature T = 1. There are three cases to distinguish, see Figure [2j 
If ab > 1/4 then the paths do not reach the hard edge at the origin. If a6 < 1/4 then 
they hit the hard edge and stick to it during a certain time interval [ti,t2]- Finally, for 
the critical separation ab = 1/4 the paths are touching to the hard edge at one particular 
time t = t* e [0, 1], given by 

t* = (1.5) 

These three cases are illustrated in Figure [5]^ a)-(c). 

In this paper we find it convenient to take fixed endpoints a,b satisfying ab = 1/4. 
We consider the time t G [0, 1] and temperature T > to be varying. Figure [3] shows 
the phase diagram in the (t, r)-plane if a = b = 1/2. The diagram shows three regions. 
Case I corresponds to temperature T < 1. Then the paths do not reach the hard edge 
and at each fixed time t G (0, 1) they are asymptotically distributed on the interval 
\p,q] = [p{t),Q{t)] with [la Sec. 4.3] 

VP :=(1 - + tVb - ^/2t{l - t)T > 0, (1.6) 
^ :=(1 - t) + tVb + V2t(l - t)T. (1.7) 

The limiting density of the paths on [p, q] is the Marcenko-Pastur density [12[ Remark 1] . 

In Cases II and III we have T > 1, with the paths reaching the hard edge in Case III 
and staying away from it in Case II. The limiting density is a transformation of the one 
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Figure 3: Phase diagram in the (t, T) plane for the non-intersecting squared Bessel paths 
with a = b = 1/2. The critical tacnode point is located at (t, T) = (1/2, 1). 

in [7] ; it is different from the Marcenko-Pastur density. We note that the limiting density 
and support are always independent of a. 

The critical curves in the phase diagram are exactly like in the corresponding model 
of non-intersecting Brownian motions [11, Sec. 1.9], thanks to [12, Remark 1]. But the 
nature of the phase transitions is different due to the presence of the hard edge. 

At the transition of Cases I and II we expect a phase transition in terms of the 
inhomogeneous Painleve II equation (|2.10p with = a + 1/2. This transition may not 
be felt in the correlation kernel but should certainly manifest itself in the critical limits 
of the recurrence coefficients of the multiple orthogonal polynomials, as in [11]. At the 
transition of Cases II and III we expect the same phase transition as in [32]. We note 
that the transition curve is given by the equation [11^ Sec. 1.9] 

a(l-t)^ + fet2 

t{l-t) ■ ^^-^^ 

The topic of the present paper is the multi-critical tacnode point {t,T) = (t* , 1) which 
lies at the transition of each of the Cases I, II and III in the phase diagram. 



2 Statement of results 

2.1 A new 4x4 Riemann-Hilbert problem 

We introduce a new 4x4 RH problem, which is a variant of the RH problem in [131114]. 
The RH problem generalizes the one in \13\ [H] in the sense that it contains non-trivial 
'Stokes multipliers' e'^'^' and e~'^^^. As in [13] we also have an extra parameter r. 

The RH problem will have jumps on a contour in the complex plane consisting of 10 
rays emanating from the origin. More precisely, let us fix two numbers ipi , cp2 such that 

0<ipi<ip2< 7r/2. (2.1) 
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Figure 4: The figure shows the jump contours Tk in the complex ("-plane and the corre- 
sponding jump matrix on T^, A; = 0, ... ,9, in the RH problem for M = M{C,). We 
denote by 0^ the region between the rays Vk and Tk+i- 



Then we define the half-lines Tk, A; = 0, . . . , 9, by 



To 

and 



(2.2) 



r5+fc = -rfc, k = o,...,4. (2.3) 

All rays Tk, k = 0,...,9, are oriented towards infinity, as shown in Figure |H We also 
denote by the region in C which lies between the rays Tk and r^+i, for /c = 0, . . . , 9, 
where we identify Fio := Fq. 

Using this definition of the rays F^, we now consider the following RH problem. 
RH problem 2.1. We look for a 4x 4 matrix valued function M{() (which also depends 
parametrically on u > —1/2 and on the complex parameters ri,r2,s,T £ C) satisfying 

(1) M(C) is analytic /or C G C \ (uLo Tfc) • 

(2) For C G Ffc, the limiting values 



M+(C)= ^ lim M(z), M_(C) 

z^Q, z on +-side oj i j. 



lim M{z) 

z—><^,z on — side o/ 



exist, where the -\--side and —-side of Tk are the sides which lie on the left and 
right ofTk, respectively, when traversing Tk according to its orientation. These 
limiting values satisfy the jump relation 



M+(C) = M_(C)Jfc(C), A: = 0,...,9, 
where the jump matrix JkiC) for each ray Tk is shown in Figure\^ 



(2.4) 
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(3) As Q ^ CO we have 



M(C) 



Ml M2 ^ / 1 



diag((-c)-^/^r^/^(-c)^/^c^/^) 

where the coefficient matrices Mi,M2, ■ ■ ■ are independent of Q, and with 



(2.5) 



A:-- 
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^2(0 = 



(2.6) 



(2.7) 



Here we use the principal branches of the fractional powers, i.e., := |C|"e*'*^''^'' 
with argument function \ argCl < vr, for any o E M. iVote that and (— C)'* /ia?;e 
a branch cut for Q along the negative and positive real line respectively. 

(4) As C ^ we have 



and 



M(c) = o(n, M-Ho = oic), 



M(c)diag(r^c^c^r'^) = o(l), 
M(c)diag(c^c-^r^n = o(l), 
Mic) = o{c-n, 



C G ill U ^^8, 
elsewhere. 



z/i/<0, 



ifu>0, 



(2. 



(2.9) 



where the 0-symbol is defined entrywise. 
More detailed information on the behavior of M(C) for (" — )• is given in Proposi- 
tion (531 

It follows from standard arguments (e.g. [8]) that the solution to the RH problem 12. II 
is unique if it exists. The existence of the solution is non-trivial and is stated in the 
following theorem, whose proof will be given in Section [6l 

Theorem 2.2. (Existence). Assume v > —1/2, ri = r2 =■ r > 0, s,t £ M. Then the 
RH problem \2. 1\ for M{C,) is uniquely solvable. The same is true when ri,r2,s,T belong 
to a sufficiently small complex neighborhood of the above values. 

The next theorem provides a connection with the inhomogeneous Painleve U equation 



2" {x) = xq{x) + 2q^{x) — 



(2.10) 



where the prime denotes the derivative with respect to x. The Hasting s-McLeod solution 
[2T] is the special solution q{x) of (j2.10p which is real for real x and satisfies 



q{x) 



q{x) 



X 
-X 



X 



X 



+00, 



-00. 



(2.11) 
(2.12) 
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Ml 



(2.15) 



We also define the Hamiltonian u{x) by 

u{x) := iq'ix)f - xq^{x) - /(x) + 2vq{x) (2.13) 

and note that 

vl^x) = -q^{x). (2.14) 

Theorem 2.3. (Mi vs. the Painleve II equation). Let the parameters v > —1/2, ri = 
r2 = 1, s,T G M in (j2.5p ~ (|2.7p be fixed. The residue matrix Mi in (j2.5p takes the form 

/a b ic id\ 

—b —a id ic 

ie if —g h 

\if ie -h gj 

where a,b,c, . . . are real valued constants depending on u, s and r. We have 

d = 2-1/3^ (2^/3(25 - r^)) , (2.16) 

c = -2-i/\(22/3(2s - T^)) + s'^, (2.17) 

with q the Hasting s-McLeod solution to the inhomogeneous Painleve II equation ()2.10p - 
()2.12p . and with u the Hamiltonian in ()2.13p . 

Theorem 12.31 win be proved in Section [5j In the homogeneous case v = Q the theorem 
recovers some results from [13| . see also [I^. In the proof we also obtain identities 
relating the other entries of Mi in terms of the entries c and d, see e.g. (|6.2p . We trust 
that the notations a, b in (j2.15p will not lead to confusion with the endpoints of the 
squared Bessel paths. 

The most difficult part of Theorem 12.31 is proving that q{x) is precisely the Hastings- 
McLeod solution to the Painleve II equation. This will require asymptotic results that 
we derive in Sections [3] and |H 

Finally, we transform the RH matrix M{C,) into a new matrix M{C,) as follows 

M(c):=diag(c^/^c-^/^c^/^c-^/^)diag((J ^'-'^^ 

The transformed matrix M(C) depends on the same parameters ri,r2, s,t,u as M(C). 
The matrix M((^) satisfies a RH problem by itself but we will not state it here. 

2.2 Critical asymptotics of squared Bessel paths 

Now we describe the critical asymptotics of the non-intersecting squared Bessel paths. 
We will work under the triple scaling limit where the endpoints a, b are fixed such that 

ab = l/4. (2.19) 

We will let the time t and the temperature T depend on n such that 

t = -^ + Kn-y\ (2.20) 
T = l + Ln"2/^ (2.21) 
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where K, L are arbitrary real constants. 

Denote again by a the parameter of the non-intersecting squared Bessel paths. It 
turns out that under the scahng hmit ()2.19p - ()2.2ip . the analysis of the squared Bessel 
paths for large n leads to the RH problem 12. H with the parameters ri = r2 = 1 and 
with V, s = s* and r = r* given by 

i/ = a + 1/2, (2.22) 

= (2.23) 

r* = -K{^+Vbf. (2.24) 

For u,v > 0, we define the tacnode kernel K'^°''^"'°'^'^{u,v; s* ,t* ,a) as 

K'^^''°'^^{u,v;s*,T*,a) = (-1 1 O) M7Hv)M+{u) (l 1 O)'^, 

27n[u — V) 

(2.25) 

where the superscript stands for the transpose. Here M((^) is defined in ()2.18p where 
the RH matrix M{() is defined with respect to the parameters ri = r2 = 1, s = s* , 
T = T* and := a + 1/2. 

The squared Bessel paths at each fixed time t G (0, 1) are a determinantal point 
process with correlation kernel Kn{x,y), see Section [2.31 below. Now we state our main 
result. 

Theorem 2.4. (Asymptotics of the correlation kernel). Consider n non-intersecting 
squared Bessel paths on the time interval [0, 1] with transition probability density (|l.ip - 
()1.4p with given starting point a > and given endpoint b > 0. Assume that ()2.19p - 
(j2.2ip hold. Then the correlation kernel Kn for the positions of the paths at time t has 
the following scaling limit as n oo with n even, 

lim -^jj^Kn (^^, = n°/2^-/2i^"^(n,^;;s*,r*,a), (2.26) 

n-^oo n^n^/^ \K'^nrl'^ K^n^l'^J 

for any fixed n, f > 0, where s* and t* are given in (|2.23p - (|2.24p and with 

K = 2(Va + \/6). (2.27) 

We expect the conclusion of Theorem 12.41 to remain valid if n — )• oo with n odd. 

Note that the factor ()2.26p has no infiuence on the correlation func- 

tions of the determinantal point process induced by i^*'^cnode^ Moreover, when compar- 
ing i^^'^cnode ^^le kernel in [13* Eq. (2.48)], note that the latter has a typo: it has 
'M-^{u)M+{vy instead of 'M-^(t>)M+(u)'. 

Remark 2.5. (Varying endpoints). Fix a*,b* > such that a*b* = 1/4 and assume that 
the endpoints a, b vary with n as 

a = a*{l + 2Lin-^/^), (2.28) 
6 = 6*(1 + 2L2n"2/3)^ (^2.29) 
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for certain real constants Li,L2. Assume again that the time t and temperature T 
vary with n as in ()2.20p - ()2.2ip . but with a,b replaced by a*,b*. Then the conclusion of 
Theorem 12.41 remains valid, with now (j2.23p replaced by 

= KHV-a + Vb)^-L + L,+L,_ ^^3^^ 

Remark 2.6. (Bessel process). Theorem 1 2 . 41 was formulated for the squared Bessel process 
(jl.ip - (jl.4p . By taking square roots, one obtains a similar statement for the (ordinary) 
Bessel process. The correlation kernel Kn{x,y) for the positions of n non-intersecting 
Bessel processes on [0, 1], starting at ^/a and ending at Vb, is expressed in terms of the 
above Kn by 

Kn{x,y) = 2y/xyKn{x^,y^), 
see |3H Remark 2.10]. The analogue of (|2.26p in this case is 

hm ^TJ^Kn f^,^) = 2u^v-''^K"''^^°''%u\v^-s\T*,a). (2.31) 

Remark 2.7. (Chiral 2-matrix model). The RH problem 1 2. II for M(Q was recently used 
to study a critical phenomenon in the chiral 2-matrix model jlOj . This leads to a critical 
correlation kernel which is essentially different from the kernel K^'^'^^°'^^ in ()2.25p and 
which is a hard edge analogue of the Duits-Geudens kernel [l3j . 



2.3 About the proof 

The positions of the non-intersecting squared Bessel paths at each particular time t £ 
(0, 1) constitute a determinantal point process. It is a multiple orthogonal polynomial 
(MOP) ensemble with respect to the weight functions ^12j 



(2.32) 



a _nx ( 2nJax\ I ^ -V^ j ( 2nJax 

wi^i{x) = x2e ^'^(^^^j, wi^2{x) = x 2 e Ttl^j^^\^—^ 

W2,i{x) = x-te-'^I, (l^) , W2,2{x) = X- Ve-n?^/,_i 



Here T > is the temperature. 

This MOP ensemble is related to the following RH problem |12j . 
RH problem 2.8. We look for a 4 x 4 matrix valued function Y satisfying 

(1) Y is analytic znC\M-|-. 

(2) For x > 0, Y possesses continuous boundary values 1+ (from the upper half plane) 
and y_ (from the lower half plane), which satisfy 

Y,ix) = Y^ix)C^ ^^^)), (2.33) 

where 1^ denotes the identity matrix of size k and W{x) is the rank-one matrix 
W{x) = (wi,i(x) wi^2{x)) (t(^2,i(a:;) ^^2,2(2;)) • 
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(3) 



As z ^ oo, z G C \ we have 



Y{z) = {I + 0(1/ z)) diag(z"i , ^ ^■ 



z 



(4) 



where ni := [n/2] and 722 := n — ni. 

y(z) has the following behavior near the origin: 



Y{z)dmg{l,l,h-\z),h-\z)) = 0{l), y-^(z)diag(/i-i(z), 1,1,1) = 0(1) 



as 2; —7- 0, z G C \ M+, where the O -symbol is defined entrywise, where the super- 
script denotes the inverse transpose and with 



There is a unique solution Y to the above RH problem. The matrix Y is constructed 
using multiple orthogonal polynomials of mixed type with respect to the modified Bessel 
weights ()2.32p : see [6, 12j. Due to the singularity of the weight matrix near the origin, 
we need the condition (4) to ensure the uniqueness of the solution. 

The correlation kernel admits the following representation in terms of Y [Uj I12j: 



Outline of the paper 

The remainder of this paper is organized as follows. In Sections [3] and [5] we calculate 
the large s asymptotics of the model RH problem for M. Section [5] proves Theorem 12.31 
on the relation of M[C,) with the inhomogeneous Painleve II equation. Section [6] proves 
Theorem 12.21 on the solvability of the model RH problem. In the final Section [7] we 
analyze the non-intersecting squared Bessel paths near the hard-edge tacnode and we 
prove Theorem 12.41 

3 Asymptotics of M{C,) for s +00 

In this and the next section we will analyze the model RH problem 12.11 for M{C,) if 
ri = r2 = 1 and r = in the limit s — )• oo. We will prove the solvability of the RH 
problem for s G M with \ s\ sufficiently large, and also establish the large s asymptotics for 
the quantities c and d in ()2.15p . More precisely, we will prove the following proposition. 
Proposition 3.1. Let v > —1/2 be fixed and suppose that ri = r2 = 1 and r = 0. Then 
for s G M with \s\ large enough, the RH problem for M(-) is uniquely solvable. 




if -1 < a < 
if a = 0, 
ifa>0. 



(2.34) 
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—i (Mi)^ ^ in (|2.15|) have the large s 

s-^+oo, (3.1) 

s +00, (3.2) 
s — )• —00. (3.3) 

Proposition 13.11 will be needed in the proofs of Theorems 12.21 and 12. 3i Note that in 
[13j we had a stronger variant of ()3.ip with 1^ = 0. This stronger variant allowed us 
to conclude directly that q in (j2.16|) is the Hasting s-McLeod solution to Painleve II. In 
contrast, the condition ()3.ip is insufficient to characterize the Hastings-McLeod solution, 
for any parameter v > —1/2, which is why we also need the asymptotics (j3.3p for 
s — )• —00. Actually, it is possible to characterize the Hastings-McLeod solution from its 
behavior for s — )• +00 alone (without the need for s — )• —00 asymptotics) |17l Sec. 11.7]; 
but the latter requires an extremely detailed asymptotic analysis that is hard to perform 
in our setting. 

The proof of Proposition 13.11 is based on the Deift-Zhou steepest descent analysis of 
the RH problem for M{C,). In this section we will perform the analysis for s — )• +cxd, in 
the next section we will consider the case where s — t- —00. 

We now turn to the analysis for s — )• +00. We will apply a series of transformations 
Ml— 7>i?i— 7>Ci— T-D, so that the matrix- valued function D tends uniformly to the 
identity matrix as s — >• +c«. The transformations will be very similar to \i2>\ Sec. 3] 
and therefore we only give a brief description. The main difference with [IS] is in the 
construction of a local parametrix at the origin with the help of modified Bessel functions; 
see Section [331 

3.1 First transformation: M ^ A 

The first three transformations will be almost identical as in [13] Sec. 3]. The first 
transformation M 1— )• ^ is a rescaling of the RH problem for M: 

9 

A(C) = dmg{s^/\ s^/\ s-^/\ s-'/')M{sO, C G C \ U r,. (3.4) 

fc=0 

The matrix A satisfies a RH problem with exactly the same jumps and behavior near 
the origin as M{(). The large C asymptotics in (|2.5p changes to 

^(c) = + y + ^ (^) ) diag((-c)-^/^ c-^/^ (-c)^/^ e^') 

X ytdiag (e-^^'^(«, e-^^'^(^), e^^'i(^), e^^"^(^)) (3.5) 

with A as in (j2.6p and with 

A = 9^iO = li-Cf^' + 2(-C)l/^ HO = IC'/' + 2C^/^ (3.6) 



Moreover, the numbers d = —i (Mi)-^ ^ and c = 
asymptotics 

d = ^ + o(.-/^), 

c = s^ + 0{s-^), 
d = ^/^ + 0(s~^), 
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Figure 5: Jump contour T,b for the RH problem for B. Note the reversion of the 
orientation of the rays Tq and (— oo,0). 

3.2 Second transformation: A^-^ B 

In the second transformation we apply contour deformations. The four rays Tk, k = 
1,4,6,9, emanating from the origin are replaced by their parallel lines emanating from 
some special points on the real line. More precisely, we replace Fi and Fg by their 
parallel rays Fi and Fg emanating from the point 2, and replace F4 and Fg by their 
parallel rays F4 and Fg emanating from the point —2. See Figure [5j 

Denote by Ej^k the 4x4 elementary matrix with entry 1 at the (j, /c)th position and 
all other entries equal to zero. We define 



A{C){I + Es^i), for C between Fi, [0, 2] and fi, 

A{C){I - £^3,1), for C between Fg, [0, 2] and fg, 

A{C){I + £^4,2), for C between F4, [-2,0] and f4, . (3.7) 

^(C)(/ - ^4,2), for C between Fg, [-2,0] and fg, 

A{(), elsewhere. 



The large C asymptotics of -B(C) are the same as for A{Q, see (|3.5p . The new jumps 
of B on the interval (—2, 2) are given by 

s+(C) = 5-(C)(/ + £^i,3), Cg(o,2), 

S+(C)=i?-(C)(/ + i^2,4), CG(-2,0). 

The jumps of B on the remaining contours are the same as those for A, provided that 
Tj is replaced by Tj, j £ {1,4,6,9}, and that the change of the orientation of f4, fg 
and (—00,0) is taken into account. (The reversion of a contour implies that the jump 
matrix is replaced by its inverse). The jump contour T,b for -B(C) is shown in Figure [5l 

3.3 Third transformation: B C 

The third transformation is a normalization of the RH problem at infinity. Define the 
'(/-functions' 

51(0 = ^(2-0'/', CeC\[2,oo), (3.8) 
52(0 = ^(2 + 0'/', C e C \ (-00, -2], (3.9) 
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with the principal branch of the power 3/2. We define the transformation B C as 

C(C) = (/ + iA(^3,i -^4,2))5(C)diag(e^^^(^),e^^'2K),e-^^i(^),e-^^2('^)), (3.10) 

where again A := 5^/2. As in [13, Sec. 3.4], one shows that C satisfies the following RH 
problem. 

RH problem 3.2. We look for a 4 x 4 matrix valued function C such that 

(1) C(C) is analytic for ( G C \ . 

(2) C has the following jumps on S^.- 

c+(c) = c_(C)Jc(C), 

where Jc is given by 



MC 
JciC 
JciC 
JciC 
JciC 
JciC 
JciC 
JciC 
JciC 
JciC 

(3) As C ^ oo, we have 

ciC) 



^2,2+-E4,4 + Sl,3-^3,l, /0rCe(2,Oo), 

(/ + e-2Ac;i{C)^^3), /orCG(0,2), 

A(9i-92)(C)^i^^VR. . - E2,i)), for C G Ta, 



il 
il 



3,4 

eA(5i-92)(C)g---i(^3^^_^2,l)), 



for C G Tg, 
for C G Ta, 



(/ + e-^^3l-92m^u^^(^E^^^ _ ^^^2))^ for C G Tr 



(/ + e-2A^;2{C)^2,4), /orCG(-2,0), 
(/ + e2A92(C)^^_2)^ /or CGf4Uf6, 

-Ell + -E'3,3 + -^2,4 - -^4,2, /or CG (-00, 



-2). 



where 



Cl = (/+iA(^3,l-E4,2)) 



diag((-c)-^/^r^/^(-c)^/^c^/^M 



Ai(I-iA(S3,l-^4,2)) + 



a! 
2 



6 




V 





a! 
2 



^A^-2^A 
6 



(3.11) 
(3.12) 
(3.13) 
(3.14) 



(3.15) 



-i\ 



\ 

-i\ 



_Ai 
2 







a! 

(3.16) 



with Ai the residue matrix in (j3.5p . 
(^J T/ie behavior of C{C) for (" — t- is the same as for M(C). More precisely, 

c(c) = o(n, c-i(c) = o(n, c^o, ifu<o, 
c(c)diag(c-^c^c^c-'^) = o(l), c^o, z/z^>o, 

where in the second line we assume that C — t- to the right o/Fa and Fg. 
As in |13] we obtain the following expressions for the numbers c and d in (|2.15p : 



(3.17) 



«\/s(Ci)i,3 + s^, d = -iVs(Ci)i,4. 



(3.18) 
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3.4 Global and local parametrices 
Global parametrix 

Away from the points 2, —2 and 0, we will approximate C by the 'global parametrix' 
C(-)(C) = diag ((2 - C)"^/^ (2 + C)-'/^ (2 - Cf/\ (2 + Cf/") A, (3.19) 

with A as in ()2.6p . and where, as usual, we take the principal branches of (2 + C)"'^''^ and 
(2 — C,Y^^ with cuts along (-co, —2] and [2,oo), respectively. See [I3i Sec. 3.6]. 



Local parametrices at the points —2 and 2 

In small neighborhoods of the points —2 and 2 we construct local parametrices C^~^^, 
C^^) to the RH problem, respectively. These parametrices are constructed with the help 
of Airy functions, exactly as in \1?>\ Sec. 3.7]. 



Local parametrix at the origin 

The next step, which has no analogue in jl3[ Sec. 3], is the construction of a local 
parametrix C'-''^ near the origin. The construction will rely on the one in (12| Sec. 5.6]. 
The local parametrix C^^^ will be defined in the disk Bp of radius p around the origin 
in C, with p > fixed and sufficiently small. It will solve the following RH problem. 
RH problem 3.3. We look for a 4 x 4 matrix valued function C^^'' such that 

(1) C7(0)(C) is analytic for ( € Bp\^B- 

(2) For C G n (Ts U Tg U Ty U Tg), we have cf^ {() = (OJciC) with Jq given in 
([XTTjl ^ dmi) . For C G (0,/9) we have 

.(0),.x _^(0),.^ / / + e-2^fiK)£;i,3, i/i/-l/2^Z>o, 



and for ^ G (— /0, 0) we have 

,(0),., „(0),., / ; + e-2*9=«)E2,4, ./,/-l/2sZZ>o, 



cr(C) = c-(Ox|;;^ (3,21) 

(^5^ Uniformly on the circle |C| = p we have for A — t- +oo i/iai 

C(°)(C) = C(°°)(C)(/ + 0(A-^)). (3.22) 

(^^J The behavior of C^^\C) for ( ^ is the same as for C{(), see (|3.17p . 
(5) C^^^ satisfies the symmetry relation 
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Figure 6: Jump contour and enclosed sectors ujq, . . . for the local parametrix C^^\ 



The idea behind ()3.20p - ()3.2ip is that we can get the correct jumps Jc as long as 
1/ — 1/2 Z>o. If the latter fails then we approximate Jc by the identity matrix for 
^ G {—p,p), but we will see that the resulting approximation error can be controlled. 

We will see below that the solution matrix C^^^ indeed exists. For now, we stress 
the symmetry relation ()3.23p and note that we have the same type of relation between 
the jump matrices on opposite rays. Then similarly as in Proposition 15.31 in Section O 
we can find the detailed local behavior of C^^^ near the origin. The outcome is that 
if 1/ — 1/2 then there exist an analytic matrix valued function E and constant 

matrices Aj. such that, for the sectors uJk in Figure El we have 



c(o)(c) = i?(c)diag(c^c^r^^'^)^fcdiag(e 

as £ ojk, k = 0, . . . ,5, with in particular 



Afli (C) ^ gAga (0 ^ g- A91 (C) ^ g- A92 (C) ^ 



(3.24) 



An 



/ 0-1 2cosz^7r\ 

2cosz^7r 1 e-"'^' 

10 

V 1 / 



(3.25) 



Note that ()3.25p is the same as the matrix Ai in (j5.7p . up to an irrelevant diagonal factor 
D. In fact the calculations in both cases are virtually the same. 

On the other hand, if — 1/2 G Z>o, then there exist an analytic matrix valued 
function E and constant matrices A^ such that 



CW(C) = S(C)i^(C)vlfcdiag(e— ,e— ,e — ,e 



A5i{C) pAgaCC) p-A3i(C) p-A92(C)^ 

and 



as C £ ujf^, k = 0, . . . ,5, where now K{() is defined in 

^0 = + + -£-3,2 + -E'4,3- 

In all cases, we have for certain constants * that (see Remark 15. 7p 



m 



(3.26) 



(3.27) 



A 





-1 






(* 


* 







1 





1 







* 


* 











1 





1 










* 


* 


Vo 


-1 





V 




^0 





* 





(3.28) 
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Now we turn to the construction of C^^\ It will be convenient to consider the following 
'square root version' dq^ of C(0): 

ci?(c) = diag(c^/^c-^/^c^/^r^/^)^diag ((| , (_\ |)) c^ic^/^) 

X diag(J, 1, 1) diag(-i, 1, 1, -i), (3.29) 

where as usual we take the principal branches of all the powers, and with J as in (|3.23|) . 
Similarly we define by replacing C(°) by C(°°) in (IMD . recall (13^9]) . We also set 

a := - 1/2, A+ := {^^ | z G Ta}, A" := {z"^ \ z £ Tg}, (3.30) 

with and A~ oriented towards the origin. Then Csq solves the following RH problem. 
RH problem 3.4. We look for a 4 x 4 matrix valued function Csq such that 



'sq [ 
.(0) 



(1) Cs7(C) is analytic for C G \ (M U A+ U A" 



(2) Ciq' has the following jumps on Bp2 r\ {M.[J 1!^+ \J i:^-): 

(Cio))+(C) = (dq°^)-(C)Jsq(C), (3.31) 

where Jgq is given by 

[ /, ifa£ Z>o, 

Jsq(C) = (/ - e^(^^-»2)K''")e-'^^(i^3,4 + i?i,2)), for C G A+ n 5^2, 
J,q(C) = (/ - e"(^^-^^)(^''')e-"-^(i?3,4 + i?i,2)), /or C G A" n 
J.q(C)=diag((_°^ J)), /orCG(V,0). (3.32) 

(3) d?(C) &e/iaues near i/ie origin — t- as follows: 

CSHC) = 0{C^% (Cg))-i(C) = 0(r/^), a<0, (3.33) 
dq°)(C) = 0(log(|C|)), (Ci°))-i(C) = 0(log(|C|)), a = 0, (3.34) 

cg)(c)diag(r/^^"/^r/^r"/') = o(i), a > o, (3.35) 

where in ()3.35p i(;e assume that C — t- m i/ie region to the right of A+ and A^ . 
(^x^j Uniformly on the circle |C| = /O^ we /lave as A — )• +00 that 

Q?HC) = Ci-)(C)(/ + 0(A-i)). (3.36) 

Proof. The jumps of dq"* on A+, A^ and easily follow from a = v — 1/2, ()3.29p and 
the jumps 

ofC7(0). Now we check the jump on M_. With the negative real axis oriented 
from left to right as usual, we have 

diag(c^/^c-^/^c^/^c-^/^)diag((J 1))^^°^^^'^'^ 



diag(c^/^c-^/^c^/^c-^/^)diag 



1 -IWI 1 

1 1 M-1 1 



idiag(-J,J) [C(°)(C'/2)]+. 
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On account of the symmetry (|3.2c{p , we obtain 



i diag {—J, J) 



By using this in ()3.29p . it is then straightforward to obtain the jump ()3.3ip - ()3.32p . 

The behavior of ciq\c,) as C — t- fohows from a = v — 1/2, (|3.29p and a careful 
inspection of (I3:2ill - (l3:28| ). □ 

The properties of ciq^ and ci^^ are exactly like those for 'T' and W^' in IJT) Sec. 5.6], 
respectively, where we are dealing with 'Case I'. We use the construction in that paper 
to build the matrix Csq solving the RH problem 13.41 see also [lOl Sec. 5.6]. 

Note that for a E Z>o the RH problem for ciq^ basically decouples into two RH 
problems of size 2x2 which can be solved using modified Bessel functions. On the other 
hand, if a Z>o then we have a genuine 4x4 RH problem and its solution requires some 
nontrivial modifications [12\ Sec. 5.6]. In the latter paper we use these modifications only 
if a < 0; but the same construction works as long as sin(7ra) ^ 0, i.e., a Z. 

Finally, we lift ciq^ back to the original setting by inverting ()3.29p : we put 

c^'\o = ^diag ((_\ I) , (\ -^)) diag(r^/^c^/^r^/^c^/^)Q(?(c^) 

X diag(i, 1, 1, i) diag( J, 1, 1), (3.37) 

if Re C > 0, and 

C7(°)(C) = diag(J,-J)C(°)(-C)diag(J,-J), ifReC<0, 

recall (j3.23p . We claim that this matrix satisfies the RH problem [33] for C^^\ It is easily 
seen that it has the right jumps. It also satisfies ()3.22p and ()3.23p . It remains to check 
the behavior of C^"-* near the origin. Recall that C(°) {() should have the same behavior 
for C — )• as C(C) in ()3.17p . Let us check this if i/ < 0. It is clear from ()3.37p that 
CW(C) = 0{C-^). Then we can again derive the formulas (j3.24p - (j3.28p . with now each 
term or (~'^ replaced by (^'^~^ or respectively. Using these formulas together 

with (|3.29p and (|3.33p we see that all the entries of E{() have a zero at (" = and so each 
term (^'^^^ or C,^'^^^ can be replaced by C,^ or C,^'^ respectively, which yields the desired 
conclusion. Similar arguments apply if > 0. 



3.5 Fourth transformation: C ^ D 



We define the fourth transformation 

r c(c)(c(-2))-i(c), 
_ c(c)(c(2))-HC), 
''^'^^ - c(c)(c(o))-HC), 
I c(c)(cM)-i(C), 



if IC + 21 <p, 
if IC-21 <p, 

if ICI < P, 
elsewhere. 



(3.38) 



Then D satisfies the following RH problem. 
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Figure 7: Contour S/) for the RH problem for D in the case — 1/2 Z>o. If u — 1/2 G 
Z>o then there is an additional jump contour on the line segment {—p,p). 



RH problem 3.5. We look for a 4 x 4 matrix valued function D such that 

(1) D is analytic in C \ S/j, where Ti£, is shown in Figure^ 

(2) D has jumps -D+(C) = D^{C)Jd{0 for C G ^d, where 

-1 

-1 
, -1 

-1 



Jd 



on IC + 2| 
on Id = p, 
on IC - 2| = p, 
on the rest ofTiQ. 



(3.39) 



(3) As C, ^ CO, we have 



(3.40) 



(4) D{Q is bounded in a neighborhood of Q = {). 
Note that, if — 1/2 Z>o, then D has no jumps in the disk around the origin. From 
the asymptotics of C{C,) and C'^'^^C,) at the origin we find that -D(C) = 0(1) as C — ^ 
(if > 0) or D{C) = OiC"^") as C (if < 0). In both cases we conclude that the 
singularity at the origin is removable and so D[C,) is analytic at = 0. 

On the other hand, if z^ — 1/2 G Z>05 then D has a jump on {—p,p) and so it is 
not analytic in the disk around the origin. But then we conclude as in the previous 
paragraph that -D(C) = 0(1) as C — ^" and so -D(C) is bounded near the origin. 

The jump matrix for D satisfies 



Jo{() = I + 0{X~^), as A ^ +00, 



(3.41) 



uniformly for ( on the circles |C| = p, \C + 2\ = p and |C — 2| = p, and the jumps on the 
remaining contours of S/j are uniformly bounded and exponentially converging to the 
identity matrix as A — )• +oo. (This holds in particular for the jump matrix on (— p, p) if 
u — 1/2 G Z>o.) By standard arguments [8l[9] we then conclude that 



DiC) = I + 



1 



as A 



,A(IC| + 1) 

+00, uniformly for ( G C \ S/j. Moreover, the matrix Di in (|3.40p satisfies 



(3.42) 



Di = KX-^ +0{X- 



X 



oo, 



(3.43) 
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for a certain constant matrix K. 

It follows from (j3.18p and the above constructions that 

c = -i^/^pi)i,3 + d=-i^s{Di)i^i, (3.44) 

where we used (j3.38p for large C, and the fact that the large C, expansion of (7*^°^^ in 
()3.19p has a diagonal form. Consequently, using ()3.43p and A = s^/"^ we get for a certain 
constant k G M, 

c = s^ + Ois-^), d= - + 0(s-^/^), s^+oo. (3.45) 

s 

3.6 Proof of Proposition 13.11 

From the above results we obtain the solvability of the RH problem for M(C) for s >> 
sufficiently large. We also obtain ()3.2p . 

Eq. (|2.16p shows that d = (2^/'^s) with q = q{s) a certain solution to the 

Painleve II equation. The asymptotics of d = d{s) for s — )• +oo in ()3.45p transform into 
similar asymptotics for q{s). Substituting these asymptotics in the Painleve II equation 
()2.10p . we find that the constant k in (|3.45p must necessarily equal iy/4. For example, if 
K > Z//4 then (j2.10p would imply that q"{s) is positive and bounded away from zero for 
all s € M large enough. This would imply that q{s) — t- cxd as s — )■ oo, which contradicts 
the fact that q{s) = 0(l/s). Similarly one can rule out the case where n < u/4. Hence 
we have k = v/A and we get ()3.ip . 

Finally, the proof of ()3.3p will be given in the next section. 

4 Asymptotics of M{() for s ^ — oo 

In this section we analyze the RH problem for M(C) if ri = r2 = 1 and r = in the 
limit where s — )• — oo, thereby completing the proof of Proposition 13.11 We will apply a 
series of transformations Mh^Ah^Bh^C>-^D, so that the matrix- valued function 
D uniformly tends to the identity matrix as s — >• — c«. The analysis will be markedly 
different from Section [3l This holds in particular for the contour deformation, the 
definition of the g-functions and the construction of the global and local parametrices. 

4.1 First transformation: M A 

The first transformation is again a rescaling of the RH problem for M. Define 

9 

A{0 = diag((-s)l/^ {-s)'/\ (-s)-l/^ (-,)-V4)M(-<), C G C \ U Ffc, (4.1) 

k=0 

where we assume that s < 0. Then A satisfies 

RH problem 4.1. We look for a 4 x 4 matrix valued function A satisfying 
(1) A{C,) is analytic for ( £C\ [jk=o^k- 
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Figure 8: Contour T,g for the RH problem for B. Note the reversion of the orientation 
of some of the rays. 



(2) A has the same jump matrix on Tj. as M . 

(3) As Q ^ CO, we have 



m = + f + ^ (^) ) diag((-c)-^/^ c"\ {-o"\ e") 



x^diag(e-^^i(^),e-^^^(«,e^^^K),e^^^K)), 



(4.2) 



with 



A = {-sfi\ uo = - 2(-c)^/^ ^2(0 = - ic'^. (4.3) 



4.2 Second transformation: B 

In the second transformation we again apply contour deformations, although in a differ- 
ent way than in Section [3j The four rays F^, A; = 2,3,7,8, emanating from the origin 
are replaced by their parallel lines emanating from some special points on the imaginary 
axis. More precisely, we replace r2 and by their parallel rays r2 and emanating 
from the point 2i, and we replace Fy and Fg by their parallel rays Fy and Fg emanating 
from the point —2i. See Figure [HI 



(J^) The behavior of A near the origin is the same as for M. 



The number d in (j2.15p now satisfies: 

d = -iv^(Ai)i,4. 



(4.4) 



We define 



B{0 



= < 



A{C,)J2 ^, for Q between F2, [0,2i] and T2-, 

1(C) J3, for C between F3, [0, 2i] and f 3, 

1(C) Jy^S for C between F7, [-2i,0] and fy, 

A{C)J^, for C between Fg, [— 2z,0] and Fg, 



(4.5) 



A(C), elsewhere, 
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Now B is analytic in C\S^, where is the contour shown in Figure [51 Note that in 

this figure we reverse the orientation on some of the rays. Then B satisfies the following 
RH problem. 

RH problem 4.2. We look for a 4 x 4 matrix valued function B satisfying 

(1) B{() is analytic for C G C \ 'Eg. 

(2) B has the following jumps on Eg; 



i?+(C) = i3-(C)^B(C), 



where is defined by 



Jb(C) = £^2,2 + ^4,4 + ^: 



1,3 



3,1, 



/orCG (0,+oo), 



JBiC) = {I + E3,i), /or CGFiUPg 
Jb(.0 = (.I + E4,2), for C^UUFe 
«^s(C) = -5-1,1 + -E'3,3 + -E'2,4 — E4^2, 



jRiC) = diag 



1 e""' 



Jj^(C) = (/ + e^^X^3,4-i^2,i)), 
^5(C) = (/- 

«^b(C) = diag ^ ^ _^ 
J5(C) = (/ + e^^X^4,3-i^i,2)), 

JgiO = (I 

(3) As C ^ 00, we have 



e-'^"^(i^4,3-i?i,2)), 



1 e""" 



forC G (-00,0) 
e"''' 

for C G fa, 

forCets, 
e-""" 

for C G fr, 
forCets. 



/or C G (0,2i), 



/or C G (-2i,0), 



m = + y + ^ (^) ) diag((-c)-^/^ c-^/^ (-c)^/^ c^/^) 

x^diag(e-^^"^(^),e-^^''^(^),e^^^^(f),e^^"^(«), (4.6) 



where A = (— s)^/^ and ^i(C); ^2(C) ^'^c given in (|4.3p . 
Now (14.41) transforms as follows: 



s(Si)i,4 



(4.7) 



4.3 Third transformation: B ^ C 

In this transformation we normalize the RH problem at infinity by means of '^r- functions'. 
The construction of the (^-functions will be markedly different from Section [3l 

Using the principal branches of the square root, we first define two functions ii{0 
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and ^2(0 by 



6(C) 



^-2-74Tc2, 6(C) 



^-2 + 74TC^, 



if Re C > 



(4.8) 



6(C) 



^-2 + V4 + C2, 6(C) 



^-2- v/ITc^, 



if Re C < 0. 



(4.9) 



Note that 




Z — 7- 00 



00. 



(4.10) 



The 5- functions are defined as the following anti-derivatives of the ^-functions: 



where 0^ is the origin reached from the first quadrant of the complex plane and from 
the second quadrant, and where the integration path of gi (or (72) is not allowed to cross 
M+ U [-2i, 2i] (or M. U [-2i, 2i] respectively). We have 



Observe that there is no integration constant in (j4.12p . Indeed by taking the limits 
z — +00 and z — t- —00 along the real line we see that the integration constant must be 
simultaneously real and purely imaginary and therefore it is zero. 

We need the following relations: 

• 9i,±(z) = 92,t(.^), for z G (-2i,2i), 

• gi,+ {x) +51 _(x) = 0, for x £ R+, 

• g2,+ {x) + g2-{x) = 0, for x G M_, 

and we also need the following inequalities for the real parts of the (7-functions: 

• Re {gi{z)) < 0, for z G Ti U Tg \ {0}, 

• Re (52(2)) < 0, for z G r4 U Fe \ {0}, 

• Re(5i,+(2) -5i,_(z)) < 0, for z G (-2i,2i), 

• Re(52,+(2) -52,-(z)) > 0, for z G (-2i,2i), 

• Re {gi{z) - g2{z)) < 0, for z G f2 U f g \ {0}, 

• Re igi{z) - g2{z)) > 0, for z G fg U \ {0}. 

Each of these relations can be straightforwardly checked from the definition of the func- 
tions 51,52- 

Now we define 

C(C) = (/-iA(S3,i -E4,2))^(C)diag(e^5i(^),e^^^2(^\e-^^^i(^), 6-^32(0). (4.13) 
Then C satisfies the following RH problem. 

RH problem 4.3. We look for a 4 x 4 matrix valued function C satisfying 




(4.11) 



5l(z) = |(-z)3/2 - 2(-z)V2 _ (-z)-l/2 + 0((_^)-3/2) 



Z — )• OO. 



z — cx) 



(4.12) 
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(1) C(C) is analytic for G C \ 

(2) C has the following jumps on T,g: 



c+(c) = C7_(c)Jg(C), 



where Jq is defined by 



E2,2 + -£^4,4 + Ei^3 - 

El^l + -^3,3 + -^2,4 - 
/gA{gi,+ (C)-gi,_{C)) 



^3,1, /or C G (0,+oo), 
/or C G TiUTg, 

forCer^UTe, 
^4,2, /or CG (-00,0) 



V 















e 

—g-^^'T* e~'^(^2,+{C)-92,-{C)) 

(/ + e'^-V'^^^^^^^-f^K))(i?3,4 - E2,i)), for C G f 2, 



















for C G r3 

\ 





_g!/7ri g-A(g2,+ {C)-ff2,-{C)) 

(/ + e''^'e-^^3iiO-92iO)(^E^^^ _ i^^^^)), /^r C G f 7, 
(/ - e-''"*e^(^'i(^)-»2(^))(^;3,4 - ^2,1)), for C G f g. 



/or C e (0,2i), 



/or C G (-2i,0), 



C CO; we have 

c{c) = + f + o (^) ) diag((-c)-^/^ c-^/^ (-c)^/^ c^/^)^. 



Now ()4.7p transforms as follows: 

d 



-i\/^(Ci)i,4- 



(4.14) 



4.4 Global and local parametrices 
Global parametrix 

The global parametrix is obtained from the RH problem for C by ignoring all the expo- 
nentially decaying entries in the jump matrices: 

RH problem 4.4. We look for a 4 x 4 matrix valued function C^°^^ satisfying 
(1) C(°°)(C) is analytic for ( G C \ (M U [-2i, 2i]). 
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(2) C(°°)(C) has the jumps 

C^^\C) = C[°^\C){E2,2 + ii^M + i?i,3 - E3,i), for C > 0, 
)(C) = ai°")(C)(i?i,i + Es,3 + i?2,4 - i^4,2), for C < 0, 

cf)(C) = ai°°)(C)diagff 'TVf ° . for Ce{0,2^), 



_^u7Ti Q y ' \^_e-'"^» 



(4.15) 



('^J As C — )• oo, we have 

a(°^)(c) = (^/+^+o(l)jdiag((-c)-^/^r^/^(-c)^/^c^/V• (4.16) 

(4) We have 

a(-)(C) = 0(|C T 2i|-V4), as C ^ ±2i, 

C(°^)(C)diag(l,C,l,C) = 0(1), as C ^ TOi/i ReC > 0. ^' ' 

Lemma 4.5. There exists a solution C^°°^ to the above RH problem. Moreover, the 
(1,4) entry of the matrix cj°°'* in ()4.16p is given by 

{C^^\^ = i. (4.18) 

Proof. We will give an explicit construction of C^°°\ Define the functions 



„.M = ^'W + f»^. = ^'<') + f?^, (4.19) 



where we recall the definition of ^1,^2 in (I4.8p - (l4.9p . Note that 

v^{z) = (-z)-i/2 + 0((-z)-5/2)^ ^^^^ (4_2o) 

z — >• 00. (4-21) 

Also define 

V?.{z) = -Vl{z), Vi{z) = -V2{z). 

We construct a compact Riemann surface IZ with four sheets T^j, j = 1, ... ,4. The 
sheets are defined by Tli = 7^3 = C\ (M+ U [-2i, 2i]) and 7^2 = 7^4 = C\ (M_ U [-2i, 2i]). 
We glue the sheets IZi and 7^3 in the usual crosswise way along the cut M+. Similarly 
we glue 7^2 and 7^4 along ]R_, and we glue TZi and 7^2 (and also 7^.3 and 7^.4) along the 
cut (— 2i,2i). We add points at infinity and at the finite branch points —2i and 2i to 
make IZ a compact Riemann surface. 

We consider the function Vk{z) to be living on the kih. sheet IZ^ of the Riemann 
surface, A; = 1, . . . , 4. Together these functions provide a bijective holomorphic mapping 
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of the Riemann surface TZ into the Riemann sphere C. The images of the four sheets of 
the Riemann surface under this mapping are given by 

viiUi) = {z G C I |z| < 1 and Rez > 0}, 

V2i'R-2) = G C I |z| > 1 and Rez > 0}, 

^;3(7^3) = {z G C I |z| < 1 and Rez < 0}, 

^4(7^4) = {z gC \ \z\ > 1 and Rez < 0}. 

The cuts of the Riemann surface are mapped into iM U T with T the unit circle. The 
images of the branch points are 

vi{±2i) = V2{±2i) = e^'^*/^, V3{±2i) = V4{±2i) = -e±'^*/^ 

wi(oo) = ^3(00) = 0, ^2(00) = ^4(00) = 00. 

The points at the origin of the four sheets also play a special role; they are mapped to 
the points ±1 and ±1 in the u-plane. 
Define the polynomial 

P{v) = 2v{v^ + 1), (4.22) 

and its square root 

defined as an analytic function in the f-plane, with cuts along the union of arcs 
iR+ U {e'^ I 9 G [-^/4, 7r/4] U [3^/4, 5^/4]}, 



and also satisfying y^P{v) > as f — )• 00 along the positive real axis. 
For = 0, we define the global parametrix C^'^\C;u) = 0(°°)(C;0) by 

C(°°)(C; 0) = {Mvk{C)))l,^, , for C G C \ Sg, (4.23) 
where the functions fj{v), j = 1, . . . , 4 are defined by 

h{v) = -^l= Uv) = ^= f3{v) = ^l= h{v) = ^=. (4.24) 
^/P{v) ^/P{v) ^/P{v) ^/P{v) 



with the above discussed branch of P{v). 

For general v, we define the global parametrix C^°°-'(C; i^) by 

C(°°)(C; v) = /CC(°^)(C; 0) diag (^e'^(^^(^», e'^^^^K))^ g/^C^aCC))^ e/^(^4(C))^ (4.25) 
where /C is a matrix of the form 

/C = / + a^3,i + /3^4,2 (4.26) 

for suitable constants a, /3. The function h in ()4.25p serves to create the jump entries 
^um fT^'^^ in (|4.15|) . Hence h must be analytic in C\T and with the counterclockwise 
orientation of T we must have 

hjf-{v) — h-{v) = ±i'7ri, for v G T±, 
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where 

T+ := Tn {t; I Imv > 0}, T_ := T n | Imv < 0}. 
We construct the function h expHcitly in a Cauchy integral form: 

h{v) = -- [ dx + - f dx. (4.27) 

2jj^v-x 2jj_v-x ^ ' 

The Plemelj formula guarantees that h has the desired jumps on T_|_ and T_. Moreover 
h tends to zero for both t; — )■ and v — )• oo. In fact, the dominant term of e'^'-^^ as 
— 7- oo is of the order 0(C~^^^) and it is then routine to find the constant matrix /C in 
gSSll-dESD so that ()iT6D holds. 

Prom the construction ()4.25p - ()4.26p we see that the (1,4) entry of c[°^^ in ()4.16p is 
independent of u. In fact, it is not too hard to calculate this entry explicitly, leading to 
the value i in (|4.18p . Finally, (|4.17p is also clearly satisfied. □ 

Local parametrices at —2i and 2i 

Near the points —2i and 2i one can again construct local parametrices C^~^*\ using 
Airy functions. We omit the details. 

Local parametrix at the origin 

At the origin we construct a local parametrix C^^^ in a similar way as in Section 13.41 
We require C^^^ to have the same jumps as C on the contour Bp n (Fi U r4 U U Fg). 
For C G (0, pi) we impose the jump 

{o)^AN ^{o)/.x / ^g(C), if - 1/2 Z>o, 



+ I 6-^^(^3,4 - E2,i) + e-^^'{Ei,2 - ^4,3), if - 1/2 G Z>o, 

(4.28) 

and for C, G {—pi,0) we put 

rW(r\ - r^^hn v / -^5(0, if - 1/2 ^ z>o, 

u+ o_ X I ^-uni^Es,^ _ E^^^^ ^ e-^\E^^2 - E,,^), ii u - 1/2 G Z>o. 

(4.29) 

We again use transformations of the type (|3.29p to move to a 'square root version' of 
the RH problems for C^^^ and C^°°\ The local parametrix in the square root setting 
can be constructed as in [12, Sec. 5.6], see also [10', Sec. 5.6], and it can then be lifted 
back to the original setting by using ()3.37p . This is very similar to Section 13.41 The 
main difference is that we are now dealing with 'Case IIF rather than 'Case F in the 
terminology of [12l Sec. 5.6]; see also Fig. [3l The details are irrelevant for our purposes 
and are omitted. 
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Figure 9: Contour for the RH problem for D in the case — 1/2 Z>o. If u — 1/2 G 
Z>o then there is an additional jump contour on the vertical line segment {—pi, pi). 



4.5 Fourth transformation: C ^ D 

We define the final transformation 



f C(C)(C(2^))-HC), if|C-2i|<p, 

a(c)(a(-2*))-Hc), ific + 2ii<p, 

C(C)(C(o))-i(C), if|CI<P, 

, C(C)(C(°^))-i(C), elsewhere. 



Then satisfies the following RH problem. 
RH problem 4.6. We look for a 4 x 4 matrix valued function D satisfying 

(1) D is analytic in C\ T,q, where is shown in Figure\^ 

(2) D has jumps D+{Q = 5_(C) ^^(C) for C G S^, where 



J 



D 



(7(0) ^(7(00 



1 

-1 
-1 



on Id = p, 
on IC — 2i\ = p, 
on IC + 2i| = /9, 
on the rest of S^. 



("5^ j4s C — >• 00, we /laue 



(J^) D(Q is bounded in a neighborhood of Q = Q. 
The jump matrix for D satisfies 



(4.30) 



(4.31) 



(4.32) 



J^(C) = / + 0(1/A), 



as A — )■ +00, 
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uniformly for z on the circles \z\ = \z + 2i\ = p and \z — 2i\ = p, and the jumps on 
the remaining contours of T,g are uniformly bounded and converging exponentially fast 
to the identity matrix as A — )• oo. We then again have estimates like in ()3.42p ~ (|3.43p . 
It follows from (|4.14p . A = (— s)^/^ and the above constructions that 

d{s) = -i^/^(cj°°^)i,4 + 0(s-^) = y/^ + 0{s-^), (4.33) 
where we used ()4.18p . This ends the proof of (|3.3p and thereby of Proposition 13.11 □ 



5 Proof of Theorem 12.3 



In this section we prove Theorem 12.31 on the connection of the residue matrix Mi with 
the Hastings-McLeod solution to Painleve II. We proceed in several steps. 

5.1 Symmetries of the Riemann-Hilbert problem for M{() 

First we will check that the residue matrix Mi can be written as in ()2.15p . To this end 
we use the symmetry relations of the model RH problem; some of these relations will 
also be used in the next sections. In what follows we use the elementary permutation 
matrix 

(5.1) 



.1 Oy 

and we use / for the 2x2 identity matrix. 

Lemma 5.1. (Symmetries). For any fixed v > —1/2 and ri,r2,s,T, we have the sym- 
metry relations 

M(C;n,r^,s,T)= (^^ ^^^M{C;ri,r2,s,T)(^^ , (5.2) 

M{-C;ri,r2,s,T)=(^^ M(C; ra, n, s, r) (^jj , (5.3) 

M-^(C;ri,r2,s,r)= (^_°^ ^) M(C; n, ra, -r) (^J "^^^ , (5.4) 
where the bar denotes the complex conjugation. 

Proof. One checks that the left and right hand sides of (|5.2p satisfy the same RH problem. 
Then ()5.2p follows from the uniqueness of the solution to this RH problem. The same 
argument applies to (j5.3p and (j5.4p . □ 

Corollary 5.2. (Symmetries of the residue matrix Mi). For any fixed v > —1/2, 
ri = r2 =: r, s,t £ M, we have 
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Consequently, Mi takes the form (|2.15p where a, b, c, d, e, f, g, h are real valued constants 
depending parametrically on v,r, s,t . If t = then we have additionally 

and then in (|2.15|) we have g = a and h = b. 



5.2 Behavior of M{() near the origin 



For the proof of Theorem 12.31 we wiU need detailed asymptotics of near the origin. 
This is provided in the following proposition. 

Proposition 5.3. (Behavior of M near the origin). Let M{Q be the solution to the 
model RH problem \2.1l Then we have, with all the branches being principal: 

• If ly— 1/2 ^ ^>o, then there exist an analytic matrix valued function E and constant 
matrices such that, for the regions in Fig. ^ 



M(c) = i?(c)diag(r,c,r",r")^fc, c&^k, k = o,...,9. 

Letting denote the jump matrix for M onV^, we have 

Afc+i = ^fcJfc+i, A; = 0, . . . ,3, 5, . . . ,9, 
where we identify Aiq = Aq, Jiq = Jq and with, among others. 



A, 



D 



/ 0-1 2cosz^7r\ 

2cosz^7r 1 e-"""' 

10 

V 1 / 



A, 



1 e 
1 
\0 1 



1 

VTvi 



-2 COS VTT 






(5.5) 
(5.6) 
1 \ 



1 / 

(5.7) 

with D = diag(e-^''*, e'"''^', -e'''^\ e^'^*). 

• If V — 1/2 ^ Z>0; then M has logarithmic behavior at the origin: There exist an 
analytic matrix-valued function E and constant matrices Ak such that 



with 



A, 



M(C) = E{OK{(:)A, 



KiC) :- 



(^i 1 1=F^ 

1 ±i 

i ±l + i 

V i / 



CGJ^fc, /c = 0,...,9, 



(C 





iClogC 







































^-u j 




(15:61) 


and 


we have. 


among others 



(5.8) 



(5.9) 



A. 



/ 1 1 \ 

10 

=Fl ±l-i ±l-i 

V ±1 ±1 / 



(5.10) 



in the case where 



zbi respectively. 
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Remark 5.4. The factorization (|5.5|) is not unique. Indeed, for any non-singular 4x4 
block diagonal matrix 

/* * 0\ 
* * 
* * 
\0 * *J 

one can substitute -E'(C) and in ()5.5p by E{C,)1C and }C~^Ak, respectively. 



K 



Proof. The proof uses ideas from [TTj, see also [221 123] among others, but the details 
will be more involved since we are now dealing with 4x4 matrices. 

First we consider the case where u — ^ Z>o- Then the matrices Ak in ()5.6p - ()5.7p 
are invertible. We can then define E by ()5.5p . i.e., we put 



E{c) = M{c)A,' diag(c-^ c~^ r , n, c g 



(5.11) 



By construction, E is analytic in C \ Ufc=o-'^fc- show that E is indeed entire. 

The relations ()5.6p show that E is analytic also on Ufc=o \ Ts- Moreover, on = M_ 
oriented from left to right we have 

E-\()E^ (c) = diag(r , r , c-^ cn+A,j,A^' diag(c-^ c-^ r , r)- (5.12) 

Recall the matrix J in (|5.ip . By a straightforward calculation. 



^4 ^5^5 ^ 



^5 {JQJ7J?,J9JoJlJ2J3JiJb) ^5 ^ 

(diag(J, -J)JiJ2hJ'iJb diag(J, -J)JiJ2-h'h-h) ^5 ^ 
(A5 diag(J, - J) Ji J2 J3 J4J5^5 ' 



(5.13) 



where in the first step we used the defining relation ()5.6p for ^fc, and in the second step 
we used that 

Jfc = diag( J, -J)Jk+5 diag( J, - J), A; = 0, . . . , 4. 



(5.14) 



To further evaluate (j5.13p . one calculates 



J 
-J 



Ji J2 J3 J4J5 



/- 














1 















-1 




-1 


V 


-1 










(5.15) 



It is easy to see that this matrix has eigenvalues e~'^^^ and —e'^^^, both with mul- 
tiplicity two. Moreover, the matrix A^ in ()5.7p was chosen in such a way that its 
first, second, third and fourth row are left eigenvectors corresponding to the eigenvalues 
e-"'^', -e"'^', -e"'^', respectively. Using this in (I5T3D we get 



(diag(e-'^"\ e"'^"*, 



diag(e-2^'^\ e-^'^^S e^'^^S e^^^^^) 

diag(C-^ C, C, C)+ diag(r , C, C , C) 
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Inserting this result in (|5.12p we see that E is analytic also on M_, and therefore in 
C\{0}. 

It remains to show that the singularity at is removable. If < then we see 
from the definition (jS.lip of -E'(C) and from the behavior (j2.8p of M near the origin that 
-E(C) = 0{C'^'^) as C — ^ 0, so (since 2v > —1) the isolated singularity at is indeed 
removable. 

If > and ^ — )• in ili then we find in a similar way (using ()2.9p and (jS.lip ) that 



i?(c) = c(c)diag(c^r^r^c 



* * * 

0** 
0** 

y* * * * 



diag(c-^ r",c,n, cej^i, 



where C(C) = 0(1) for ^ — )• and we used the particular form of the matrix A^^, 
cf. ()5.7p . It follows that -E'(C) is bounded near the origin in Qi and so cannot be a pole. 
Since cannot be an essential singularity either, the singularity is indeed removable. 
This ends the proof of the proposition if u — ^ ^ Z>o- 

Next we consider the case where u — ^ € Z>o- In this case, ()5.1ip and ()5.12p are 
replaced by 

E{C) = M{C)A^^K{C)-\ C^^k, 



and 



for C G 



Now ()5.15p reduces to 



(5.16) 
(5.17) 



-J 



/t« 

1 

-1 =Fi 

V-1 



\ 

±i 

-1 
^2iJ 



(5.18) 



in the case where 
same value e~^'^^ = - 



= ibz respectively. All the eigenvalues of this matrix have the 
'* = =pz. The eigenvector space corresponding to this eigenvalue 
is two-dimensional. The matrix A^ in ()5.10p was chosen in such a way that its third and 
first row, and similarly its fourth and second row, form Jordan chains in the sense that 



A. 



-J 



'h'h'h'h'hA^ ^ 



/i 1 0\ 
i 1 
i 



Hence by (|5T3]l . 















yo 












1 


^\ 


2 


/-I 





2i 


o\ 





i 





1 







-1 





2i 








i 













-1 





Vo 












Vo 








-1/ 
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A straightforward calculation shows that this matrix can be written in the form 

recall (|5.9|) . Inserting this in (jS.lZp shows that E is analytic on By definition, E is 
also analytic in C \ M_ . Finally, one shows as before that the singularity at the origin is 
removable, so E is analytic in the entire complex plane. □ 

We need the following symmetry relation for -E'(C)- 
Lemma 5.5. (Symmetry of E). Assume that ri = r2 =: r and define J in (15. ip . With 
the notations in Proposition 15.31 we have the symmetry relation 

^(0= (o _°j)^(-C)diag(l, 1,-1,-1), CgC. (5.19) 

Proof. We will give the proof for — ^ Z>oj the other case follows from similar 
considerations. 

Prom the definition ()5.5p of E{C,) and the symmetry relation ()5.3p of the RH problem 
for M(C) it follows that 

E{o = diag( J, -j)E{-c) diag((-c)^ (-c)^ (-c)-^ {-cm 

X ^fcdiag(J,-J)^^_^5diag(^^c-^c^c), (5.20) 

if C £ for any = 0, . . . , 4. We will use this relation for A; = 0. Then 

^0 diag( J, -J)A^^ = ^5 {JqJjJsJqJo diag( J, - J)) A^^ 
= A5 (diag(J, -J) Ji J2 J3^4 J5) ^5 ^ 
= diag(e-^^^e-'^^^-e'^^^-e^^^), (5.21) 

where we used again the relations ()5.6p . ()5.14p and the fact that the first, second, third 
and fourth row of A^ are left eigenvectors corresponding to the eigenvalues e~'^^^, e~'^^^, 
^_^vm q£ matrix ()5.15p . respectively. Inserting ()5.2ip in ()5.20p (with /c = 0), 
and using the principal branches of C and C,~'^ with C, £ $^5, we obtain the desired 
relation (|5J9]l for C ^ ^5- From (f5H]l and (fOO]) it fohows that the relation (f5J9]l 
remains valid in the other sectors k = 0, ... ,4 and hence in the entire complex 

plane. □ 

We are now ready to prove the following lemma. 
Lemma 5.6. Assume that ri = r2 =: r. With the notations in Proposition \5.3l we have 

— M-HC) = ^^"'(0 + diag(l, 1, -1, -1)E-\C), C G C, (5.22) 

if u — ^ 1i>o- The first term in the right hand side of (|5.22p is analytic, whereas the 
second term behaves as 

7diagff° ''^,(' -:))+Oil), C-O. (5.23) 



If u — ^ G Z>o then the same conclusions hold except that, in the right hand side of 
(j5.22p . we have an extra term of order 0{z'^'^~^) as z ^ 0. 
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Proof. We restrict ourselves to the case where v — ^ ^ ^>oj the other case is similar. 
The equation ()5.22p follows immediately from ()5.5p . Next, we will show that the second 
term of (j5.22p behaves as (j5.23p near C = 0- To this end, we must show that 

i?(0)diag(l,l,-l,-l)i?-i(0)= ((J (5.24) 

This follows by substituting C = in dEED- □ 

Remark 5.7. Substituting = in ()5.19p we get a number of equations for the entries 
of E{Q). It is then straightforward to check that -E(O) must be of the form ()3.28p . 



5.3 Lax pair and compatibility conditions 



Now we proceed with the proof of Theorem l2.3l We will derive the Lax pair correspond- 
ing to the RH problem for M(C). Throughout this section we take a fixed u > —1/2, 
T S M and we let ri = r2 = 1. First we obtain a differential equation with respect to C. 
Proposition 5.8. We have the differential equation 

dM 



UM 



where ( recall ()2.15p ) 



U :-- 



( -C + T 

-d+vC^ 
-i{-C + g + a + s) 
\ -i{h + h) 



d + vC^ 

C — T 
-i(b+h) 

-iiC + g + a + s) 



Proof. On account of ()2.5p . we have for ^ — )• oo, 





C + T 

-d-uC'^ 





i 

d-uC^ 
—c — tJ 



(5.25) 



(5.26) 



dM 



M- 



^ Ml M2 








i(l-<-i) 





\ 





— r 





i{l + sC 




i{C-s) 





r 







V 


-i{C + s) 





— T 


/ 



M2 



+ ... +0(r^). (5.27) 



Ml Mf - 

Since the jump matrices in the RH problem for M are independent of C, the matrix 
^^^^M{Q~^ is analytic for ^ G C \ {0}, and it has a first order pole at ( = given 
by (j5.23p . The terms in (j5.27p which are polynomial in (" are as in (j5.26p and so by a 
standard application of Liouville's theorem we obtain ()5.25P " (|5.26p . □ 



From the proof of Proposition 15.81 we can deduce more. As already mentioned, we 
know that the C~ ^-coefficient of ()5.27p must have the form in ()5.23p . Evaluating the top 
rightmost 2x2 block of this ^"^-coefRcient we find after some calculations. 



a + c^ - d'^ + g - 
-b + h + 2Td 

Hence we obtain the following result. 



b-h-2Td 
-c^ + d^ - g + s 



0. 
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Lemma 5.9. The numbers a, b, c, d, 



in (|2.15|) satisfy the relations 



g = —a — (? + dF + s, 
h = b- 2Td. 



(5.28) 
(5.29) 



Remark 5.10. By subtracting the (1,2) and (3,4) entries of the C, ^-coefficient in ()5.27p . 
we obtain the additional relation 



/ 



-26c + c^d + 2Tcd + 2T^d - d^ 



2sd + V. 



(5.30) 



Next we obtain a differential equation with respect to s: 
Proposition 5.11. With the notations (j2.15p . we have the differential equation 



dM 
ds 



V -.= 2 



/ c d -i 0\ 

d c i 

i(— C + a + (?) i{b — h) —c —d 

\ i{h-b) -i{C + a + g) -d -c) 



(5.31) 



Proof. The matrix ^^^^^^ M(C)~'^ is again analytic for C G C \ {0}. It is also analytic 
at C = 0, due to the fact that the matrices in (|5.5p and (|5.8p are independent of s. 
Moreover, (|2.5p yields 



dM 
ds 



■M- 



r Ml 











-2i 


0\ 













2i 




-2< 











v 





-2iC 





0/ 



Ml 



+ ... 



(5.32) 



as C — )• oo. The terms which are polynomial in C, are given by V in ()5.3ip and so by 
Liouville's theorem we obtain (j5.3ip . □ 

Now we turn to the compatibility condition of the two differential equations (j5.25p 
and (|5.3ip . By computing the partial derivative = in two different ways, 

thereby making use of the relations ()5.25p and (|5.3ip . one obtains the matrix relation 



dU dV 

— = — + VU-UV. 
OS oC 



(5.33) 



Writing this matrix relation in entrywise form leads to the following lemma. 
Lemma 5.12. The numbers b,c,d in (j2.15p . viewed as functions of s, satisfy the fol- 
lowing system of coupled first order differential equations: 



c 
d' 

y 



4d^ + 2s, 

4{-b + cd + Td), 

-46(c + r) + 4d(c2 + 2rc + 2t'^) - 4d^ - 6sd + 2i/, 



(5.34) 
(5.35) 
(5.36) 



where the prime denotes the derivative with respect to s. 
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Proof. One can write (j5.33|) in its entrywise form with the help of (|5.26|) and (j5.31|) . This 
leads to a system of 4 x 4 = 16 coupled first order differential equations. After some 
straightforward calculations, one checks that these 16 relations reduce to 4 independent 
conditions. The (1, 1) matrix entry yields (|5.34p . while the (1,2) and (3,2) entries yield 



d' = A{cd-Td-h), (5.37) 
b' + h' = -4:{2d{a + g) + bc+ih-b)T + ch + sd-i^). (5.38) 

The relations ([SSSD-dSSSD now foUow from dOTjl - ffOHD and <^M)-<^M)- □ 
Lemma 5.13. The number b in ()2.15p satisfies 

b' = {c + T)d' + Ar'^d - - 6sd + 2u, (5.39) 

where the prime denotes the derivative with respect to s. 

Proof. This follows from substituting ()5.35p in (j5.36p . □ 
Lemma 5.14. The number d in ()2.15p satisfies 

d" = 32d^ + 32sd - IGr^d - 8z^, (5.40) 

where the prime denotes the derivative with respect to s. 

Proof This follows from dOijl . (l5:35]l and (fOOD . □ 



Proof of Theorem 12.31 It is easily seen that any solution d = d{s) to ()5.40p is of 
the form 

d = 2-i/3g(^22/3(2s-r2)), 

where g is a solution to the Painleve II equation (|2.10p . This yields (j2.16p . 

Next we must show that q{s) in (|2.16p is indeed the Hasting s-McLeod solution to the 
Painleve II equation, i.e., we must establish ()2.1ip - ()2.12p . If r = then these relations 
follow from the asymptotic behavior of d in (j3.ip and (j3.3p . If r G M \ {0} one could 
perform a similar analysis as in Sections [3] and H] to show that the same asymptotics 
hold true. An alternative proof follows from Section 16.31 There we will use the Lax pair 
relations, with q the Hastings-McLeod solution, to define the matrix M{C,) if r 7^ 0, and 
we will then prove that the resulting matrix M{C,) indeed satisfies the RH problem 12. II 

Finally, we need to prove the formula ()2.17p for c. By comparing ()5.34p with ()2.14p . 
we obtain the following relation for the derivative of c = c{s) with respect to s: 

c'(s) = -2-V3^(^(22/3(2s _ r^))) + 2s. 

Then ()2.17p follows by integrating this relation with respect to s. The fact that the 
integration constant in (|2.17p is zero, follows from (j3.2p if r = 0. The case where r 7^ 
can be obtained as in the previous paragraph. □ 

As in |14] . we can also obtain a differential equation with respect to the variable r. 
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Proposition 5.15. With the notations (|2.15|) . we have the differential equation 
dM 



dr 



WM, 



W :-- 



1 c 


-2b 





-2id\ 


-2b 


-c 


2id 








-2if 


c 


-2h 


Vif 





-2h 


-c ) 



(5.41) 



The compatibility condition of the differential equations ()5.25p and ()5.4ip now reads 

dU _ dW 



+ wu - uw. 



(5.42) 



For later use, we note that the (1,2) entry of (j5.42p yields, with the help of (|5.30p . 

— ^ = b-cd-Td. (5.43) 
4r OT 



6 Proof of Theorem 12.2 



In this section we prove Theorem 12.21 on the solvability of the RH problem 12. II for M(^) 
for real values of ri = r2 =: r, s, r. The proof uses Fredholm operator theory and the 
technique of a 'vanishing lemma' if r = 0. For t ^ we will follow [14:\ Sec. 5]. 



6.1 Fredholm property 

In this section we show that for any ly > —1/2 and ri,r2,s,r G C with ri,r2 having 
positive real part, the singular integral operator associated to the RH problem 12.11 is 
Fredholm with Fredholm index zero. This follows from a powerful general procedure 
that has been applied in different settings in the literature [THl E51 ES] • For completeness 
we give a brief account of those arguments that involve the particular structure of our 
4x4 RH problem. Hereby we closely follow the exposition in [23^ Sec. 2.3]. 

Let ID) := {C G C I Id < 1}. Denote the matrix in ([23]) by 

HiC) := diag((-C)-^/^ C'^\ {-Cf'\ e'')A diag (e-^^(^)+<, e'^^^^-^^, e^^K)+^^ , e^^^^" 

In what follows we assume for convenience that v — 1/2 Z>o. We define a new 4x4 
matrix valued function M^^\C,) by 

M(i)rn = / ^(CM7'diag(C-^C-^C^n, Ce%nD, 

j = 0, . . . , 9, with Aj the matrices in (j5.5p . Here if' stands for the complement of the 
closed unit disk. 

By Proposition Owe have that M^^) (C) is analytic in D. Let S := IJ^^g U and 
orient it as in Figure [TOl Then S is a complete contour, in the sense that C \ S allows a 
decomposition as the disjoint union of two sets: C \ S = U r2_, $1+ n Vt- = 0, such 
that E is the positively oriented boundary of 17+ and the negatively oriented boundary 
of Denote = Q.j n as shown in Figure [TUl 

Now M^^\C,) satisfies the following RH problem. 
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Figure 10: Contour E = \J%o(^j U Ej) for the RH problem for M^^\(). The figure also 
shows the decomposition of C \ E into two disjoint regions and 

RH problem 6.1. We look for a 4 x 4 matrix valued function M^^\Q (which also 
depends parametrically on v > —1/2 and on ri,r2, s,t & C) satisfying 

(1) M(i)(C) is analytic forCGC\E. 

(2) (0 = M«(C)Jm(i)(C) for C e S, where 

Jmw{C) =h{OJmjH-\c), C g TjnW, j {0,5}, 

Jjv,(i)(C)=diag(C,C,r",r")vli(C)^-nC), CeS,, J even, 

JM(l)(c)=^(c)^7'(c)diag(c-^^^c^^, ces,-, j odd, 

(3) M«(C) = / + 0(C-^) for C ^ 00. 

We observe that J]^(i) — I decays exponentially as C — t- oo along S. Also observe that 
the contour S has 11 points of self- intersection, of which one is the origin and the other 
10 lie on the unit circle. At each fixed point of self-intersection lying on the unit circle, 
say P, order the contours that meet at P counterclockwise, starting from any contour 
that is oriented outwards from P. If we denote the limiting value of the jump matrices 
over the jth contour at P by Jjvf(i) j{P)^ j = 1) ■ ■ ■ ) 4, then a little calculation shows that 
we have the cyclic relation 

JmW,i{P)JmW,2^P)JmW,z{P)JmW^^{P) = I- 

The situation at the origin is trivial since all jump matrices there are the identity. 
Now we will obtain a factorization 

J^a)(C)=: (t;_(C))-V(C)- 
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Outside small neighborhoods of the 10 points of self-intersection on dD we choose the 
trivial factorization f+ = V- = I. Using the cyclic relations we are then able 

[I8l [35] to choose a factorization of J^ii) (C) in the remaining neighborhoods in such 
a way that u+ (or w_) is continuous along the boundary of each connected component 
of rj+ (or respectively). Standard arguments [181 l35] then imply the Predholm 
property with Fredholm index zero. 



6.2 Existence of M(C) if r = 

Now we show the existence of Af(^) if ri, r2 > 0, s G M and r = 0. Thanks to the 
Fredholm property in the previous section, the existence follows if we can prove that the 
'homogeneous version' of the RH problem 12. H obtained by replacing the /+0(C~^) series 
in ()2.5p by 0(C~^), has only the trivial solution M = 0. This approach is known as a 
vanishing lemma [HI [181 [35]. In the present setting it can be established in essentially the 
same way as in \1'6\ Sec. 4]. There are a few minor differences due to the z^-dependence 
of the jump matrices, and the fact that there is singular behavior at the origin. 

For example, defining the products of jump matrices 



J 



M(4) 



J8J9J0J1J2 



( 

€ 

V 



1 







1 

— V'l 

1 





i 

) 



X G 



/I e"""' - e"^'^* 






Vo 



1 







1 



1 / 



X G 



then a key property in the proof of the vanishing lemma is that, with 



-h 



G(C) :=diag(e''l^^^e■ 
J{x) :=e:i(x)J,,,(4) G+(x) 

we have the rank-one property 

J(x) + J^(x) = 2ei^(x) (1 e~ 

where the superscripts ^ and ~^ denote the conjugate transpose and the inverse con- 
jugate transpose, respectively. Using these formulas, it is straightforward to adapt the 
proof of the vanishing lemma in Sec. 4]. □ 



CeC, 

X E M, 



0) (1 e'''^' 0) 61^ (x), X G 



Remark 6.2. An alternative approach for proving the existence of M if ri = r2 = 1, 
s G M and r = 0, avoiding the vanishing lemma, is to use an analogue of the approach 
from [m Sec. 5] that is outlined in Section 16.31 below. We should then keep r = 
fixed and work with the system of differential equations (|6.ip , with the second equation 
replaced by ()5.3ip . From Proposition 13. 1 1 we already know that the RH matrix M exists 
for s G M large enough. The outlined approach then allows to extend this conclusion to 
ah s G M. 
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6.3 Existence of M(C) if r 7^ 



(6.1) 



Finally we prove the existence of M(Q if ri = r2 =: r and s, r G M with r 7^ 0. To this 
end we use the approach of Duits and Geudens [HI Sec. 5]. By rescaling we may assume 
that ri = r2 = 1. We consider the system of two differential equations 

f^(C,r) = C/(C,T)M(C,r), 
^{C,T) = W{C,T)MiC,T), 

with U and W as in (|5.26p and (j5.4ip respectively, where we now define 

d = 2-1/3^(22/3(25 _ ^2))^ 

C=-2-V3^,(22/3(2s-t2))+s2, 
b=^^ + cd + Td, 

f = 2T^d -^^-c'^d-d^-2sd + u, 
h = 4zM + Jd-Td, 



2) 



g + a = — + d"^ + s, 



where q is defined as the Hastings-McLeod solution to Painleve II and u is the corre- 
sponding Hamiltonian. The formulas in (j6.2p are well-defined for every r E M since 
the Hastings-McLeod function q has no poles on the real line [5j. These formulas are 
compatible with the ones we already established. Indeed, for the entries c and d this 
follows from Theorem 12.31 The formula for b follows from ()5.43p . and the formulas for 
f,g, h in ([62]) are then obtained from dOSjl - ffOO]) . 

We note that (j6.2p is exactly the same as ^14^ Eq. (5.5)] except for the extra term v 
in the formula for / and in (|2.1Up . 

The next two lemmas are proved in exactly the same way as in [14] . 
Lemma 6.3. (of. J14\ Lemma 5.1]). The matrices U and W defined above satisfy the 
compatibility relation ()5.42p . 

Lemma 6.4. (cf. 114[ Lemma 5.2]). Fix r G M and let U be as defined above. Let S be 
one of the four complex sectors 

Si = {C G C I < argC < 27r/3 + e}, E2 = -E^, E3 = -Ei, E4 = E^, 



where e > is fixed and small and where the bar denotes the complex conjugation. Then 
the equation ^(C^t) = U{C,t)N{(^,t) has a unique fundamental solution N in the 
sector E with the following asymptotics as (" — )• 00 within this sector, 



N{0 



^ ] ) diag((-C) 



X A diag \^e 

with A,0i,92 given in ()2.6p - (j2.7p . and with Ni of the form 



iVi 



i{C)+rC e-e2(C)-< e^ico+rc g92{C)-TC 



(6.3) 



/* 


b 


* id\ 


-b 


* 


id * 


* 


if 


* h 


Vf 


* 


-h */ 



The asymptotics in (|6.3p are uniform within any closed sector of ( away from 
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Following [14, Sec. 5] , the fundamental solutions to the differential equation ^ = UN 
in the four sectors in Lemma[63]can be used to construct a matrix function M = M{(^, r) 
satisfying the correct jumps and asymptotics for ^ — )• oo in the RH problem 12.11 and 
satisfying also the second differential equation in ()6.ip . Using the latter differential equa- 
tion, we find that the matrices E{Q in Proposition 15.31 satisfy ^ = WE. Since W is 
analytic and uniformly bounded and we already know that E{Q is analytic for r = 
(due to Proposition 15.31 and the existence established in Section [6.2p . we then find that 
E{Q is analytic for all r G M. This shows that M has the correct behavior at the origin 
in the RH problem O □ 



7 Proof of Theorem 12.4 



In this section we prove Theorem 12.41 on the critical behavior of the non-intersecting 
squared Bessel paths near the hard-edge tacnode. The proof will follow from a Deift- 
Zhou steepest descent analysis of the RH problem 12.81 for Y{z). The main technical step 
will be the construction of the local parametrix near the origin and this is where the 
model RH problem 12. 11 for M(^) will be used. 

In what follows we consider n non-intersecting squared Bessel paths with fixed end- 
points a,b satisfying ()2.19p . We assume that n is even. We study these paths at an 
n-dependent time t and temperature T which vary according to the triple scaling limit 

7.1 Modified A- functions 

First we define some auxiliary objects that will be needed during the steepest descent 
analysis. Inspired by [131 Sec. 6.3], define 

a := (1 - t) + tVb - y/2t{l - t)T, (7.1) 

/3 := (1 - t)y/^ + tVb + y/2t{l - t)T, (7.2) 

and 

7 := (a + /3 + 2^02 + ^2 _ q,/3^ /3, (7.3) 
S := (^a + P - ^/(^^TW^^) (7.4) 

A little calculation shows that under the triple scaling assumptions (|2.19p - (|2.2ip . we 
have for n — )• 00 that 

7 = ^ + 0(,.-./3), (") 
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Lemma 7.1. Sec. 6.3]). There exists an analytic function A : C \ [0, oo) — t- C such 

that 

0, xe(o,7], 



ReA+(x) = ReA_(x) 
ImA+(x) = — ImA_(2;) ■ 



> 0, xe(7,oo), 



>0, xG(e,7], 
_= 7r/2, X G (7,00), 

with e a number that goes to zero as 0(n^^/^) as n 



(7.6) 
(7.7) 



00; moreover 



^^^^ 4t(l-t)T 2tT^ 2(1 -t)r" 
for a certain constant £, and 



1 



log{-z)+i + Oiz-'), z^oo, (7.8) 



X{z) 



-X-zfl' + {-zf"G{zl 



t{l-t)T 

where G{z) is an analytic function in a neighborhood of z = which satisfies 

V7 , ^,^-2/3n 



G(0) 



3t{l - t) 



+ 0(n-^/^), 



n 



00. 



(7.9) 



(7.10) 



Proof. Take \{z) to be the function which is called Ai(2;) in [131 Sec. 6.3] with oi := ^/a, 
bi := y/b and pi := ^. Note that the temperature T = 1 in [T3] but the extension to a 
nontrivial temperature is straightforward. □ 

Now we construct the functions 



Ai(z) 

A2(^) 

A3(^) 



2A(-\/^) + 2t{l-t)T 

2A(yi) + 2t(i^j)T ~ 
-2A(Vi) + 



az 



tT + (l-t)T 



(l-t)T 



2t(l-t)T 
-2A(-^) + 2t(l-t)T 



vl^ _| v6£ 1-2/' + TT? 



(7.11) 



(l-t)T 



for iblmz > 0. 

Lemma 7.2. The \-functions have the following asymptotics for z — )• cxd.- 



A2(-) = TiT^ - ^ + 4/ - 1 log z + ^ + O(z-i), 
A3(^)= (f^ + ^log^-^+0(z-i), 
A4(^) = - (f^ + i log z + ^ + O(z-i), 



/or certain constants c\ and C3. 
Proof. Immediate from the definitions. 



□ 
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Note that the asymptotics in Lemma 17.21 are exactly the same as for the A- functions 
in [12^ Lemma 4.8]. Also note that 

^ ^ -2K{^+Vbfn-'/' + 0{n-^/'), n ^ oo, (7.12) 



tT (1 - t)T 

due to our triple scaling assumptions ()2.19p - ()2.2ip . 
In what follows we write, with 



g:=7'. (7.13) 

Lemma 7.3. We have 

'^i,±(a;) = A2,=p(x) =F vri, x G M_, 

A2,±(a:) = A3,=p(x), xe[0,g], 

^3,±ix) = A4^q:(x) ±TTi, X G M_, 

and 

Re (A2(x) — A3(x)) > 0, x e {q, oo). 

There is a Jordan curve A2 running from to q in the upper half place so that 

Re(A2(z)-A3(z)) <0, zeAf\B„ 

where A2 = and Bs denotes the disk around the origin with radius p > sufficiently 
small. Similarly, there is a Jordan curve A'^ running from oo to in the second quadrant 

of the plane so that, with A~[ = A'f and p > sufficiently small, 

Re(Ai(z)-A2(z)) >0, Re(A3(z) - A4(z)) > 0, zeAf\Bp. 

Proof. The equalities follow from the definitions. The inequalities are intrinsic in our 
definition of A in the proof of Lemma |7.H with the help of [13, Lemmas 7.8, 7.9]. □ 

An illustration of the Jordan curves and A^ in shown in Figure [TTJ 



7.2 The transformations Y t-^ X t-^ U t-^ T 

In \T2\, it is shown how to apply a steepest descent analysis to the RH problem 12.81 for 
Y{z). The authors apply there a series of transformations y 1— )• X 1— )• f7 1— t- T. We can 
apply exactly the same transformations here except that we replace the A-functions in 
|12j by the modified A-functions ()7.1ip . The resulting matrix T satisfies the following 
RH problem (see Case I in [121 Proposition 5.4]): 

RH problem 7.4. We look for a 4 x 4 matrix valued function T{z) satisfying 
(1) T{z) is analytic for z G C\ (M U A^ U A^) , where the contours are as in FiaureUli 
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Figure 11: Jump contours in the RH problem for T{z) in the analysis of the non- 
intersecting squared Bessel paths in the multi-critical case. 



(2) For each of the oriented contours in Figure [771 T has a jump T+(z) = T^{z)Jt{z) 
where 

Jt(x) = diag ( 1, ( ° /orxe(0,g), 



Jy(x) = / + e-"(^2,+ {a;)-A3,+ W)_g^^^^ xeR+\ [0, q], 

Jt{x) = diag {^(y^^ J) ' (y^^ J)) . forx£ (-oo,0), 
J^{z) =1 + e"(^2W-A3(.))^3^2, for z G , 

where the notation Ej^^ defined in Section [Y. 
(3) As z ^ oo we have 



r,.) = (/ + o(i))diag,.-v.,../.,,-v.,,./.,_L,i,,((i ;),('; 



uniformly for z G C \ M. 
(4) T{z) behaves for z — )■ as 

r(z)diag(|z|"/M2|-"/Mz|"/Mz|-"/2)) = 0(1), if a > 0, 

r(z)diag((log|z|)-i,l,(log|z|)-i,l) = 0(1), ifa = 0, (7.14) 

T{z) = 0(W2), T-\z) = 0(W2), i/ -1< a < 0, 

where we assume that z — t- in i/ie region between Af and A^. 



7.3 Global parametrix 

The global parametrix (z) is the solution to the RH problem obtained by setting all 
the exponentially small entries in the RH problem for T(z) equal to zero (see Lemma [7.3|) : 
RH problem 7.5. We look for a 4 x 4 matrix valued function P^°°\z) satisfying 

(1) P^°°\z) is analytic for zeC\ (M_ U [0, q]). 
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(2) For X £ M_ U [0, q] we have the jump 



pl-^>{x) = pr'{x)di^s(i,( \ J 1,1), 



-1 



Pl^>{x)=P'_'{x)di8.g 



(3) As z ^ oo we have uniformly for z G C \ M, 

1 



p(oo)(. 



1 
-1 



1 



X G (0,g), 
, j;GM_. 



/ + O ( ^ j j diag(^-l/^ z^/^)-^ dia^ 



1 A A i 

i 1 j ' U 1 



f^J p(°°)(z) has at most fourth-root singularities at the special points q and 0. 

The global parametrix was constructed in a more general setting in [121 Sec. 5.4]. In 
the present case the construction can be made fully explicit: 
Lemma 7.6. The solution to the RH problem \ 7. 5| for is given explicitly by 



p(-)(.) = Cdiag(z-V4,,i/4,,-i/4^,i/4)_l^diag((^ ;),(^ \ 



X ^-idiag(72-i(V^),7r'(^/^),7i(V^),72(^/i))A, (7.15) 



where 



A:-- 



1 



/I OA 



V2 



1 
-i 1 
Vi 1/ 

C:= 1+^(^2,3-^4,1), 7i(^) := f:-^V^\ I2{z) 



Proof. Straightforward verification. 



(7.16) 

1/4 

(7.17) 
□ 



7.4 Local parametrix at the point q 

Inside a fixed small disk around the point q we construct a local parametrix P^'^\z) to 
the RH problem with the help of Airy functions. Such a construction is standard and 
well-known and we do not go into the details. 



7.5 Local parametrix at the origin 

In this section we construct a local parametrix near the origin. To this end we will use 
the model RH problem for M(C). The construction will also use the 'squaring trick' of 
Its et al. [23]. 
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Figure 12: Jump matrices in the RH problem for = N{Q. 



Transformation of the Riemann-Hilbert problem for M{() 

Recall the model RH problem 12. II for Af(C). We put := a + 1/2 and we set 

iV(C):=I)-Miag((j j),l,l)M(C)diag((j J)) ' ^' ^' 



where 



D = diag(— i, 1, 1, — i). 



The jumps for are shown in Figure [T2l 

The asymptotics of as C — ^ oo is given by 



/ + f + f + o (1) ) diag(c-^/^ (-c)-^/^ (-c)^/^ c^/^) 

X Idiag fe-^^(«-<,e-^i(f)+^^,e'^i(^)+<,e^^(^)-< 



with A as in ()7.16p . The behavior around infinity can be rewritten as 



1 



mo = diag(c-^/^ (-c)-^/^ (-c)^/^c^/^)-4 [i+^ + ^+o ^^.^^^ 

X diag (e-'^(0-rC^^-eAC)+rC^^e,(0+rC^^e,iO-rC 



with the lb sign as C — ^ c« within the upper/lower half plane. Here 

Ni,± = diag (l, e^''^/^ 0, o) A^i diag (o, 0, e^''*/^ l) A. 



(7.18) 



(7.19) 



(7.20) 
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For further use, we record the symmetry relation 
N{-C;ri,r2, s,t) = diag(^(^^^ o) ' (-1 

xiV(C;r2,ri,s,r)diag((° "^^^ , "p^)) . (7.21) 

This easily follows from (j5.3|) . Note that the order of ri and r2 differs on both sides of 
the equality. 

Construction of the local parametrix 

Now we construct the local parametrix P^'^^ around the origin. In [13) we constructed 
P^'^^ inside a shrinking disk Bp of radius p = n^^^^ around the origin, see also [H]. In 
the present setting, we will be able to work inside a fixed disk Bp with radius p > fixed 
but sufficiently small. The fact that we have a fixed (rather than a shrinking) disk will 
greatly simplify some of the technical details. 

RH problem 7.7. We look for a 4 x 4 matrix valued function P^^^ such that 

(1 ) {z) is analytic for z e Bp\ {RU Af U Af). 

(2) For z G n (M U A]^ U A2 ), -P^°^ has the same jumps as T, see RH problem\T^ 

(3) Uniformly on the circle \z\ = p we have for n — )■ cxd that 

p(0)[z) = P^'^\z){I + 0(n-i/3)). (7.22) 

(4) The behavior of P^^\z) for z ^ is the same as for T{z), see ()7.14p . 

We need some functions f{z), ri(z), r2{z) and constants s, r. We will use a signifi- 
cantly simpler construction that in [13j and [14J. Note that there are many parameters 
f,ri,r2,s,T and there is a certain freedom in how to define them. 

We set (recah ^17171 ) 

K = 2(V^ + Vb), f{z) = K^z, ri(z) = 3G(V^)k-3/2, (7.23) 



r2(z) = ZG{-^z)K 



'3/2 „ _ "V ^ -1/2 



-1/2 _ ( Vb 



(7.24) 



t{l-t)T ' \tT {l-t)T 

recah ([72]). Thanks to ([ZSD and (mTTl we have 

-|r2(z)/(z)3/2 - 2./(z)i/2 _ rf{z) = -Ai(z) + h{z), 

-|ri(z)(-/(z))3/2 - 2s(-/(z))V2 + rf{z) = -^z) + h{z) ^ vri, 

|ri(z)(-/(z))3/2 + 2s{-f{z))^l^ + Tf{z) = -A3(Z) + h{z) ± TTi, 

|r2(z)/(z)3/2 + 2s/(z)i/2 _ rfiz) = -\^{z) + h{z), 

with h{z) := 2t{i-t)T + 

The following limits exist: 

:= li,„ n^r^s = li,„ „V^^f-''f = '<HV-a+f,'-L_ 

n— ^-oo >oo t[\ — t) 2 

T* := lim n^/V = -K{^/^+Vbf. 
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The first formula follows from (j2.19|) - (|2.21|) and (|7.5|) . The second one follows from 
()7.12p . From ()7.10p we also find that for z in a neighborhood of the origin, 

hm r,{z) = \; + 0(zi/2) = 1 + 0(zi/2), j = 1, 2. 

n-s-oo — t) 

Now we define the local parametrix P^^\z) at the origin by 

= ^(z)iV(n2/3/(z);ri(z),r2(z),n2/3s,n^/V)A(z), (7.26) 
with N{C) defined in frT8|) . and with 

i^(.) = P(-)(Z) diag(C^/^ (-C)'/^ (-C)-'/^ r^/^), C := n2/=^/(z), (7.27) 
A(z) = exp (-n (^-^-^-^ + 2^) ) diag(e"^^(^), e"^^(^), e"^^(^), e"^^(^)). (7.28) 

Lemma 7.8. The matrix-valued function E{z) in (j7.27p zs analytic for z G C\IR_. For 
X < /ias the jump 

i^+(x) = i^_(x)diag((_°^ J)), x<0. (7.29) 

Proof. It follows from the definition ()7.27p of E{z) and the fact that P^°°\z) is analytic 
for z G C \ M that E{z) is also analytic for z G C \ M. For x G M+, one checks that the 
jumps of and of the rightmost factor in (j7.27p cancel each other out so that E{x) 
is analytic for x G M+. A similar calculation yields the jump (j7.29p of E{x) for x G 
(Alternatively, the lemma could be proved with the help of the explicit formula ()7.15p .) 

□ 

We claim that the matrix P^^\z) in (|7.26p satisfies the RH problem 17.71 The jumps 
of P^^'^ on M-i- U Af U are easily checked. We are assuming here that the curves Ag", 
Af, A^ , Ag are chosen near the origin so that they are mapped to the rays Fi, F2, Fg, 
Fg respectively under the map z 1— t- y/z. For x G M_ we have 

Pf^\x) = ^+(x)iV+(n2/3/(a;);ri(x),r2(x),n2/3s,n^/3r)A+(x) 

= ii;-(^)((_°i j))iV+(n2/3/(a:);ri(x),r2(x),n2/3s,ni/3^)A+(x) 

= ii;_(x)Af_(n2/3/(x);ri(x),r2(x),n2/3s,nV3r)((_°^ o) ' {-I o)) ^^^'^^ 

- ^°'(^)-.((.°, ;).(_», ; 

where in the second line we used (|7.29p and in the third line we used (j7.2ip and the fact 
that ri{z) = r2{—z). This yields the required jump on M_. 

The matching condition (f?:^ fohows from dZ^H]), fTMp and fT^ . The fact that 
P'^*')(z) has the correct behavior for z — t- follows from a = v — 1/2, ()7.26p . ()7.18p and 
a careful inspection of (|5.5p - (|5.1Up and (|3.28p . 
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Figure 13: Jump contour S/j in the RH problem for R{z). 



7.6 Final transformation of the Riemann-Hilbert problem 

Using the global parametrix p(°°) and the local parametrices P^^-* and P^^\ we define 
the final transformation T i— )• by 

r r(z)(pW)-i(^), if kl <p, 

R{z)=l r(z)(pW)-i(z), ii\z-q\<p, (7.30) 
[ r(z)(p(°°))-i(z), elsewhere. 

Lemma 7.9. With the contour S/j as in Figure [T3. we have R-^- = R^Jr on S/j where, 
for a suitable constant c > 0, the jump matrices behave as 

Jr{z) = I + 0{n^^), on the circle \z — q\ = p, 
Jr{z) = I + 0{n~^/^), on the circle \z\ = p, 

/g-cni/3\ 

Jr{z) = I + O \ \ I j ' uniformly for z on the other parts of Tir. 

Proof. Only the statement about the circle |2:| = p requires further clarification. This 
follows from (|7.22p . Recall that we have here a disk with fixed radius p, in contrast to 
the shrinking disk in |131 [Hj . □ 

Note that R{z) has no jumps in the disk around the origin. From (|7.30p and the 
behavior of T{z) and P^^'> (z) for z — t- 0, we also find that R{z) can have no pole at z = 
and therefore the singularity at the origin is removable. Hence R{z) is analytic at z = 0. 

By standard arguments, see e.g. [8], we obtain from Lemma 17.91 the estimate 
as n — )• oo, uniformly for z G C \ Tjr. 



7.7 Proof of Theorem EH 

In this section we prove Theorem 12.41 We start from the expression for the kernel in 
(j2.35p . Unfolding the transformations Y ^ X ^ U i— s-Twe obtain from [12, Eq. (6.1)] 
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that 

Kn{x,y) = (0 -e"^2.-(y) ^nX.^.iy) Q)T-\y)T+{x) 

1tti[x — y) 

X (O e-"^2,+W g-nA3,+(x) Q^r^ (7 32) 

If X, y lie within radius p of the origin, then it fohows from T{z) = R{z)p(^\z) and 
the definition fL^Eh of that 

i^n(x,y) = ^^^-^(0 -1 1 0)iV+\n2/3/(y);n2/3s,nV3r)i?;i(y)i?;\y) 

X i?+(x)S+(x)iV+(n2/V(2;);n2/3s,n^/3^) (0 1 1 O)^. (7.33) 

Now we put 

for fixed u,v > 0. Using (j7.23p it then follows that 

n2/3/(x) ^ ^/^, n2/3/(y) ^ x/^, 
and we also recall from (j7.25p that 

Furthermore, we have from (|7.3ip and Cauchy's formula that 

R-\y)Rix) = 1 + (:^) = / + O (^) . (7.35) 

The constants implied by the O-symbols are independent of u and v when u, v are 
restricted to compact subsets of M4. . 

Next we estimate the matrix E{z) in ()7.27p . We claim that the transformed matrix 



^(.):=E(.)diag((^ ;),(^ ;))diag(r^/^c^/^^^/^c^/^), c 

(7.36) 

is analytic near z = 0. Indeed it has no jumps, by virtue of Lemma l7.8| and moreover 
it behaves as 0{z~^^^) as z — t- so there is no pole at z = 0. 

From dZSS]), (I73i]l and (17:2711 we get the estimates E{x) = 0{n^/^) and E-^{y) = 
0(?7-^/^). By the analyticity of E{z) and Cauchy's formula we then find 



E-^{y)E{x) =1 + (n(x - y)) = I + O 



u — V 



as n — )• 00, uniformly for u,v £ M"*". Combining this with ()7.35p we find 

hm E-\y)R-\y)R+{x)E{x)=I. 
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To use the above estimate, we should first express the matrix E in (|7.c{3|) in terms 
of its transformed counterpart E in (|7.36p . This substitution releases an extra factor 
which multiplies from the left or right the matrix N or N^^, respectively, in (j7.33p . By 
combining this with the above estimates we find 

n^oo n-^n^i-i \K^n^l-i n^n^/-^J 2m{u — v) ^ ' ^ 

X iV+(V^;s*,r*) (0 1 1 O)^, (7.37) 

with 

^ {c}!^) := diag(c^/^c-^/^c^/^c-^/^)diag , -;)) n (c^/^) , 

for C = u,v. Equivalents, by (I7T8]1 and ([238]) . 

iV (C^/') = diag(-i, 1, 1, i)M(C) diag [{^_. J) , 1, -i) . 
Inserting this in (j7.37p we get 

hm =^^^^^^(-1 1 ^)M-\v-s\t*) 

n^ca K.^n^l-i \K^r&l-^ K^n^i'^J 2m(u — v) ^ ' ^ 



X M+(n;s*,r*) (1 1 O)^. 



This finally proves (j2.26p . 
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